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Introduction 



The description of neutrino interactions with nuclei, and nuclear matter in general, 
is relevant to the study of many different problems, from supernovae explosions to 
neutron star cooling, as well as to the determination of the properties of neutrino 
itself, most notably its mass. 

The appearance of a supernova is the last stage of the evolution of stars with 
initial mass bigger than ~ 4 M©, where Mq k, 2 xlO^^ g denotes the solar mass 
[H [2|, [3] . Although the first attempts to simulate a supernova explosion date more 
than 30 years back [1], the problem is not solved yet. In fact, the results of numerical 
calculations predict that, due to loss of energy carried away by neutrinos, produced 
in the dissociation of atomic nuclei in the core, the explosion does not occur. 

The systematic uncertainty associated to simulations depends heavily on the val- 
ues of the neutrino-nucleon and neutrino-nucleus reaction rates used as inputs. As 
many existing programs use values obtained from models based on somewhat over- 
simplified nuclear dynamics [5], El [7] , such uncertainty may be significantly reduced 
adopting more realistic models, which have proved very successfuU in the description 
of electro-magnetic interaction of nuclei (see, e.g., Refs. [SEI). 

A similar problem arises in the field of neutrino physics, which is rapidly develop- 
ing after the discovery of atmospheric and solar neutrino oscillation [TOj [TTl [12] . The 
experimental results point to two very distinct mass differences 0, Am2„; ^ 8.2xl0"5 
eV^ and |Am^j^| 2.5 x 10"^ eV^. Only two out of the four parameters of the 
three-family leptonic mixing matrix Upmns are known: 612 ~ 34° and 623 ~ 45°. 
The other two parameters, ^13 and 6, are still unknown: for the mixing angle ^13 
direct searches at reactors [15] and three-family global analysis of the experimental 
data [16], [17] give the upper bound Oi^ < 11.5°, whereas for the leptonic CP- violating 
phase 6 we have no information whatsoever. Two additional discrete unknowns are 
the sign of the atmospheric mass difference and the ^23-octant. 

Neutrino oscillation experiments measure energy and emission angle of the charged 
leptons produced in neutrino-nucleus interactions, and use the obtained results to 
reconstruct the incoming neutrino energy. Hence, the quantitative understanding 

third mass difference, Am|^^p ~ 1 cV^, suggested by tlic LSND experiment [13], has not 
been confirmed yet [l4| . 
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of the neutrino-nucleus cross section, as well as of the energy spectra and angular 
distribution of the final state leptons, is critical to reduce the systematic uncertainty 
of data analysis. A number of theoretical studies aimed at providing accurate pre- 
dictions of neutrino- nucleus scattering observables are discussed in Refs. [T8 1 IT9 1 I20]. 

It is important to realize that, while neutrinos interacting in stellar matter typ- 
ically have energies of the order of few MeV, the energies involved in long baseline 
oscillations experiments are much larger. For example, K2K takes data in the region 

= 0.5-3 GeV. 

The huge difference in kinematical conditions is reflected by different reaction 
mechanisms. For neutrinos of energy < 10 MeV, since the spatial resolution of 
incoming particle is much bigger than the average distance between nucleons, the 
nuclear response is largely determined by many-body effects. 

On the other hand, at energies of the order of 1 GeV, it is reasonable to expect 
that the scattering process on a nucleus reduce to the incoherent sum of elementary 
processes involving individual nucleons. Furthermore, in this kinematical regime one 
needs to take into account the fact that elementary neutrino-nucleon interactions 
can give rise to inelastic processes, leading to the appearance of hadrons other than 
protons and neutrons. 

The results of electron- and hadron-induced nucleon knock-out experiments have 
provided overwhelming evidence of the inadequacy of the independent particle model 
to describe the full complexity of nuclear dynamics [22l [23]. While the peaks cor- 
responding to knock-out from shell model orbits can be clearly identified in the 
measured energy spectra, the corresponding strengths turn out to be consistently 
and sizably lower than expected, independent of the nuclear mass number. This 
discrepancy is mainly due to the effect of dynamical correlations induced by the 
nucleon-nucleon force, whose effect is not taken into account in the independent 
particle model. 

Nuclear many body theory provides a scheme allowing for a consistent treatment 
of neutrino-nucleus interactions at both high and low energies. Within this approach, 
nuclear dynamics is described by a phenomenological hamiltonian, whose structure 
is completely determined by the available data on two- and three-nucleon systems, 
and dynamical correlations are taken into account. 

Over the past decade, the formalism based on correlated wave functions, orig- 
inally proposed to describe quantum liquids [21], has been employed to carry out 
highly accurate calculations of the binding energies of both nuclei and nuclear mat- 
ter, using either the Monte Carlo method [251 EEl [27] or the cluster expansion for- 
malism and the Fermi Hypernetted Chain integral equations [281 1291 [30] . 

A different approach, recently proposed in Refs. [311 [32] exploits the correlated 
wave functions to construct an effective interaction suitable for use in standard 
perturbation theory. This scheme has been employed to obtain a variety of nuclear 
matter properties, including the neutrino mean free path [21] and the transport 
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coefficients [32| [33]. 

In this work we describe the apphcation of the formahsm based on correlated 
wave functions and the effective interaction to the calculation of the weak response 
of atomic nuclei and uniform nuclear matter. 

The Thesis is structured as follows. 

In Chapter [H we briefly describe the main features of nuclear matter and nuclear 
forces. 

In Chapter [21 we focus on the theory of nuclear matter, introducing the correlated 
states and the cluster expension formalism, needed to define the effective interaction. 
The applicability of perturbation theory within this framework is also discussed. 

Chapter [3] is devoted to an overview of the many body theory of the nuclear 
matter response, and its connection to the spectral function. 

In Chapter m we discuss the assumptions underlying the impulse approximation, 
and its applicability in the high energy regime. After a short description of the 
main features of the spectral function obtained from nuclear many-body theory, 
we compare the calculated electron-nucleus cross sections to data, and show the 
predictions of our approach for charged current neutrino-nucleus interactions. 

Finally, in Chapter [5] we focus on the low energy regime. We develop a formalism 
based on the effective interaction of Ref. [32] and an effective Fermi transition 
operator, obtained at the same order of the cluster expansion. The effects of both 
short and long range correlatios, described within the correlated Hartree Fock and 
Tamm-Dancoff approximations, are discussed. 

Note that in this Thesis we use a system of units in which h = h/27r = c = 1, 
where h is the Planck constant and c is the speed of light in the vacuum. 
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Chapter 1 

Nuclear matter and nuclear forces 



Nuclear matter can be thought of as a giant nucleus, with given numbers of protons 
and neutrons interacting through nuclear forces only. As the typical thermal energies 
are negligible compared to the nucleon Fermi energies, such a system can be safely 
considered to be at zero temperature. 

A quantitative understanding of the properties of nuclear matter, whose calcula- 
tion is greatly simplified by translational invariance, is needed both as an interme- 
diate step towards the description of real nuclei and for the development of realistic 
models of matter in the neutron star core. 

The large body of data on nuclear masses can be used to extract empirical infor- 
mation on the equilibrium properties of symmetric nuclear matter (SNM), consisting 
of equal numbers of protons and neutrons. 

The (positive) binding energy of nuclei of mass number A and electric charge Z, 
defined as 

= 1 [Zmp + (A - Z)mn + Zm^ - M(Z,A)] , (1.1) 

where M(Z,A) is the measured nuclear mass and rrip, m„ and me denote the proton, 
neutron and electron mass, respectively, is almost constant for A> 12 , its value 
being ~ 8.5 MeV (see Fig. II. ip . The A-dependence is well described by the semi- 
empirical mass formula 

_ = _ [ - aA''' - - + A ap-^ ] . (1.2) 

The first term in square brackets, proportional to A, is called the volume term and 
describes the bulk energy of nuclear matter. The second term, proportional to the 
nuclear radius squared, is associated with the surface energy, while the third one 
accounts for the Coulomb repulsion between Z protons uniformly distributed within 
a sphere of radius R. The fourth term, that goes under the name of symmetry energy 
is required to describe the experimental observation that stable nuclei tend to have 
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the same number of neutrons and protons. Moreover, even-even nuclei (i.e. nuclei 
having even Z and even A — Z) tend to be more stable than even-odd or odd-odd 
nuclei. This property is accounted for by the last term in the above equation, where 
A —1, and -|-1 for even-even, even-odd and odd-odd nuclei, respectively. Fig. 11.11 
shows the different contributions to i?(Z,A)/A, evaluated using Eq. (11. 2p . 

In the A oo limit, and neglecting the effect of Coulomb repulsion between 
protons, the only contribution surviving in the case Z = A/2 is the term linear in 
A. Hence, the coefficient ay can be identified with the binding energy per particle 
of SNM. 
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Figure 1.1. Upper panel: A-dependence of the binding energy per nucleon of sta- 
ble nuclei, evaluated according to Eq. (jl.2p with ay = 15.67 MeV, Og = 17.23 MeV, 
Oc = .714 MeV, ttA = 93.15 MeV and Op = 11.2 MeV. Lower panel: the solid line 
shows the magnitude of the volume contribution to the binding energy per nucleon, 
whereas the A-dependence of the surface, coulomb and symmetry contributions are 
represented by diamonds, squares and crosses, respectively. 



The equilibrium density of SNM, po, can be inferred exploiting the saturation of 
nuclear densities, i.e. the experimental observation that the central charge density 
of atomic nuclei, measured by elastic electron-nucleus scattering, does not depend 
upon A for large A. This property is illustrated in Fig. II. 2[ 

The empirical values of the binding energy and equilibrium density of SNM are 

po = O.iefm^-'^ , E = -15.7MeV. (1.3) 
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In principle, additional information can be obtained from measurements of the exci- 
tation energies of nuclear vibrational states, yielding the (in)-compressibility module 
K. However, the data analysis of these experiments is non trivial, and the resulting 
values of K range from K ~ 200 MeV (corresponding to more compressible nuclear 
matter, i.e. to a soft equation of state (EOS)) to ~ 300 MeV (corresponding to 
a stiff m^) [31]. 

The main goal of nuclear matter theory is deriving a EOS at zero temperature 
(i.e. the density dependence of the binding energy per particle E = E{p)) capable 
to explain the above data starting from the elementary nucleon-nucleon (NN) in- 
teraction. However, many important applications of nuclear matter theory require 
that its formalism be also flexible enough to describe the properties of matter at 
finite temperature. 

Unfortunately, due to the complexity of the fundamental theory of strong inter- 
actions, the quantum chromo-dynamics (QCD), an ah initio description of nuclear 
matter at finite density and zero temperature is out of reach of the present com- 
putational techniques. As a consequence, one has to rely on dynamical models in 
which nucleons and mesons play the role of effective degrees of freedom. 

In this work we adopt the approach based on nonrelativistic quantum mechanics 
and phenomenological nuclear hamiltonians, that allows for a quantitative descrip- 
tion of both the two-nucleon bound state and nucleon-nucleon scattering data. 

In this Chapter we outline the main features of nuclear interactions and briefly 
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describe the structure of the NN potential models employed in many-body calcula- 
tions. 

1.1 Nuclear forces 

The main features of the NN interaction, inferred from the analysis of nuclear sys- 
tematics, may be summarized as follows. 

• The saturation of nuclear density (see Fig. 11.21) . i.e. the fact that density 
in the interior of atomic nuclei is nearly constant and independent of the 
mass number A, tells us that nucleons cannot be packed together too tightly. 
Hence, at short distance the NN force must be repulsive. Assuming that the 
interaction can be described by a nonrelativistic potential v depending on the 
inter-particle distance, r, we can then write: 

v{r) > , |r| < rc , (1.4) 

Tc being the radius of the repulsive core. 

• The fact that the nuclear binding energy per nucleon is roughly the same for 
all nuclei with A> 12 suggests that the NN interaction has a finite range ro, 
i.e. that 

v{r) =0 , |r| > ro . (1.5) 

• The spectra of the so called mirror nuclei, i.e. pairs of nuclei having the same 
A and charges differing by one unit (implying that the number of protons in a 
nucleus is the same as the number of neutrons in its mirror companion), e.g. 
ifN (A = 15, Z = 7) and {A = 15, Z = 8), exhibit striking similarities. 
The energies of the levels with the same parity and angular momentum are 
the same up to small electromagnetic corrections, showing that protons and 
neutrons have similar nuclear interactions, i.e. that nuclear forces are charge 
symmetric. 

Charge symmetry is a manifestation of a more general property of the NN inter- 
action, called isotopic invariance. Neglecting the small mass difference, proton and 
neutron can be viewed as two states of the same particle, the nucleon (N), described 
by the Dirac equation obtained from the Lagrangian density 

C = i,N{il^d^-m)ijN (1.6) 

where 

= ( ' ) • 
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p and n being the four-spinors associated with the proton and the neutron, re- 
spectively. The lagrangian density (11. 6p is invariant under the SU (2) global phase 
transformation 



where a is a constant (i.e. independent of x) vector and the tj {j = 1,2,3) are 
Pauli matrices (whose properties are briefly collected in Appendix El)- The above 
equations show that the nucleon can be described as a doublet in isospin space. 
Proton and neutron correspond to isospin projections +1/2 and —1/2, respectively. 
Proton-proton and neutron-neutron pairs always have total isospin T=l whereas 
a proton- neutron pair may have either T = or T = 1. The two- nucleon isospin 
states |T,Mj') can be summarized as follows (see also Appendix R|) 



Isospin invariance implies that the interaction between two nucleons separated by 
a distance r = |ri — and having total spin S depends on their total isospin T 
but not on its projection M^- For example, the potential v{r) acting between two 
protons with spins coupled to S = is the same as the potential acting between a 
proton and a neutron with spins and isospins coupled to S* = and T = 1. 

1.1.1 The two-nucleon system 

The details of the NN interaction can be best understood in the two-nucleon system. 
There is only one NN bound state, the nucleus of deuterium, or deuteron (^H), 
consisting of a proton and a neutron coupled to total spin and isospin S = 1 and 
T = 0, respectively. This is a clear manifestation of the spin dependence of nuclear 
forces. 

Another important piece of information can be inferred from the observation 
that the deuteron exhibits a non vanishing electric quadrupole moment, implying 
that its charge distribution is not spherically symmetric. Hence, the NN interaction 
is non-central. 

Besides the properties of the two-nucleon bound state, the large data base of 
phase shifts measured in NN scattering experiments (the Nijmegen data base [35] 
includes ~ 4000 data points, corresponding to energies up to 350 MeV in the lab 
frame) provides valuable additional information on the nature of NN forces. 



(1.8) 



|1,0) 

1,-1) 
|0,0) 
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The theoretical description of the NN interaction was first attempted by Yukawa 
in 1935. He made the hypothesis that nucleons interact through the exchange of a 
particle, whose mass fi can be related to the interaction range tq according to 

ro ~ - . (1.9) 
Using To ~ 1 fm, the above relation yields /i ~ 200 MeV (1 fm~^ = 197.3 MeV). 




Figure 1.3. Feynman diagram describing the one-pion-exchange process between 
two nucleons. The corresponding amplitude is given by Eq. (|1.10p . 



Yukawa's idea has been successfully implemented identifying the exchanged par- 
ticle with the TT meson (or pion), discovered in 1947, whose mass is i^ttS^ 140 MeV. 
Experiments show that the pion is a spin zero pseudo-scalar particlqj (i.e. it has 
spin-parity 0~) that comes in three charge states, denoted vr"*", vr" and 7r°. Hence, it 
can be regarded as an isospin T=l triplet, the charge states being associated with 
isospin projections Mt=+ 1, and —1, respectively. 

The simplest vr-nucleon coupling compatible with the observation that nuclear 
interactions conserve parity has the pseudo-scalar form ig'j^T, where g is a coupling 
constant and r describes the isospin of the nucleon. With this choice for the in- 
teraction vertex, the amplitude of the process depicted in Fig. 11.31 can readily be 
written, using standard Feynman's diagram techniques, as 

{f\M\t) = -ig ,2_ 2 ■ ^2) , (1-10) 

^The pion spin has been deduced from the balance of the reaction tt"*" ^ p + p, while its 
intrinsic parity was determined observing the 7r~ capture from the K shell of the deuterium atom, 
leading to the appearance of two neutrons: tt^ + d ^ n + n. 
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where m,r is the pion mass, k = p[ —pi = p2 —P2, k'^ = k^k^ = k^ — |kp, u{p,s) is the 
Dirac spinor associated with a nucleon of four momentum p = {p,E) (E=a/p^ + m^) 
and spin projection s and 

(ti ■ Ta) = vl'TV2 Vv'TVi , (1-11) 

rji being the two-component Pauh spinor describing the isospin state of particle i. 

In the nonrelativistic limit, Yukawa's theory leads to define a NN interaction 
potential that can be written in coordinate space as 

-.2 ^—mwr 



£^(n-T2)(^i-V)(^2-V) 

9^ ml 1 



(47r)2 4m2 3 



(cTi ■ 0-2) + S 



12 



3 3 
- + — 

X X'^ 



e 

X 



where x = I r I and 



— [a,-cT2)5^'\v)\ , (1.12) 



3 

•^12 = ;^('^i • r)(cr2 ■ r) - (cti • 0-2) , (1.13) 



reminiscent of the operator describing the non-central interaction between two mag- 
netic dipoles, is called the tensor operator. The properties of S12 are summarized in 
Appendix El 

For g'^/{47r) ~ 14, the above potential provides an accurate description of the 
long range part (|r| > 1.5 fm) of the NN interaction, as shown by the very good fit of 
the NN scattering phase shifts in states of high angular momentum. In these states, 
due to the strong centrifugal barrier, the probability of finding the two nucleons at 
small relative distances becomes in fact negligibly small. 

At medium- and short-range other more complicated processes, involving the 
exchange of two or more pions (possibly interacting among themselves) or heavier 
particles (like the p and the u mesons, whose masses are rup = 770 MeV and m^^ = 
782 MeV, respectively), have to be taken into account. Moreover, when their relative 
distance becomes very small (|r| < 0.5 fm) nucleons, being composite and finite in 
size, are expected to overlap. In this regime, NN interactions should in principle 
be described in terms of interactions between nucleon constituents, i.e. quarks and 
gluons, as dictated by QCD. 

Phenomenological potentials describing the full NN interaction are generally 
written as 

v = v^ + vr (1.14) 

where is the one-pion-exchange potential, defined by Eqs. ( 11.12^ and ( ]1.13p . 
stripped of the 5-function contribution, whereas vr describes the interaction at 



11 



1 - Nuclear matter and nuclear forces 



medium and short range. The spin-isospin dependence and the non-central na- 
ture of the NN interactions can be properly described rewriting Eq. fll.Mp in the 
form 

^(^j) = [vTs{rij) + SsiVtTirij)Si2] P2S+in2r+i , (1-15) 

ST 

S and T being the total spin and isospin of the interacting pair, respectively. In 
the above equation P2S+1 {S = 0,1) and II2T+1 (T = 0,1) are the spin and isospin 
projection operators, whose definition and properties are given in Appendix S 

The functions VTs{.fij) and Vtririj) describe the radial dependence of the inter- 
action in the different spin-isospin channels and reduce to the corresponding com- 
ponents of the one-pion-exchange potential at large r^j. Their shapes are chosen in 
such a way as to reproduce the available NN data (deuteron binding energy, charge 
radius and quadrupole moment and the NN scattering data). 

An alternative representation of the NN potential, based on the set of six oper- 
ators (see Appendix E]) 

Of' = [l,(o-, ■ a^),Si,] ® [1,(t, ■ Tj)] , (1.16) 

is given by 

vi^j) = i2''^''\r,,)0^^ . (1.17) 

n=l 

While the static potential of Eq. fll.l7p provides a reasonable account of deuteron 
properties, in order to describe NN scattering in S and P wave, one has to include 
the two additional momentum dependent operators 

0:;=''« = L ■ S ® [1,(t. ■ T,)] , (1.18) 

L being the orbital angular momentum. 

The potentials yielding the best available fits of NN scattering data, with a 
X^/datum ~ 1, are written in terms of eighteen operators, with 

^„=9,...,i4 ^ [L2,L2(a,.^,.),(L-S)^]®[l,T,-T,] , (1.19) 

CTj,Sij]®T,j , {T,i + T,j) (1.20) 

r){Tj-v)-{n-Tj) . (1.21) 

The tak.e care of small charge symmetry breaking effects, due to the 

different masses and coupling constants of the charged and neutral pions. 



where 
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Figure 1.4. Comparison between the ^Sq and phase shifts resulting from the 
Nijmegen analysis [35] (open circles) and the predictions of the Argonne vis (AVI 
8) and one-pion-exchange (OPEP) potentials. 



The calculations discussed in this Thesis are based on a widely employed poten- 
tial model, obtained within the phenomenological approach outlined in this Section, 
generally referred to as Argonne vis potential [36]. It is written in the form 



18 



;i.22) 



n=l 



As an example of the quality of the phase shifts obtained from the Argonne 
Vis potential, in Fig. 11.41 we show the results for the ^Sq and partial waves, 
compared with the predictions of the one-pion-exchange model (OPEP). 

We have used a simplified version of the above potential, obtained including the 
operators O^-^, originally proposed in Ref.j37]. It reproduces the scalar part of the 
full interaction in all S and P waves, as well as in the ^Di wave and its coupling to 
the ^Si. 

The typical shape of the NN potential in the state of relative angular momen- 
tum i = and total spin and isospin S = and T = 1 is shown in Fig. II. 5[ The 
short-range repulsive core, to be ascribed to heavy-meson exchange or to more com- 
plicated mechanisms involving nucleon constituents, is followed by an intermediate- 
range attractive region, largely due to two-pion-exchange processes. Finally, at large 
interparticle distance the one-pion-exchange mechanism dominates. 
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Figure 1.5. Radial dependence of the NN potential describing the interaction 
between two nucleons in the state of relative angular momentum i = 0, and total 
spin and isospin S = and T = 1. 

1.1.2 Three-nucleon interactions 

The NN potential determined from the properties of the two-nucleon system can be 
used to solve the many-body nonrelativistic Schrodinger equation for A > 2. In the 
case A = 3 the problem can be still solved exactly, but the resulting ground state 
energy, Eq, turns out to be slightly different from the experimental value. For exam- 
ple, for ^He one typically finds Eq = 7.6 MeV, to be compared to E^xp = 8.48 MeV. 
In order to exactly reproduce Eexp one has to add to the nuclear hamiltonian a 
term containing three-nucleon interactions described by a potential Vijk- The most 
important process leading to three-nucleon interactions is two-pion exchange asso- 
ciated with the excitation of a nucleon resonance in the intermediate state, depicted 
in Fig. Ol 

The three-nucleon potential is usually written in the form 



where the first contribution takes into account the process of Fig. 11.61 while is 
purely phenomenological. The two parameters entering the definition of the three- 
body potential are adjusted in such a way as to reproduce the properties of and 
^He [38]. Note that the inclusion of Vijk leads to a very small change of the total 
potential energy, the ratio {vij)/{Vijk) being ~ 2 %. 

For A > 3 the Scrodinger equation is no longer exactly solvable. However, very 
accurate solutions can be obtained using stochastic techniques, such as variational 
Monte Carlo (VMC) Green function Monte Carlo (GFMC) [25]. 



V^,k = V'l + Vi 



ijk ' 



(1.23) 
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N N N 




N N N 

Figure 1.6. Diagrammatic representation of the process providing the main con- 
tribution to the three- nucleon interaction. The thick sohd hne corresponds to an 
excited state of the nucleon. 



The GFMC approach has been successfully employed to describe the ground state 
and the low lying excited states of nuclei having A up to 10. The results of VMC 
and GFMC calculations carried out using realistic nuclear hamiltonian, summarized 
in Fig. 11.71 show that the nonrelativistic approach, based on a dynamics modeled 
to reproduce the properties of two- and three-nucleon systems, has a remarkable 
predictive power. 




Figure 1.7. VMC and GFMC energies of nuclei with A < W compared to 
experiment (from Ref . [39j ) . 
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Chapter 2 

Nuclear matter theory 



Understanding the properties of matter at densities comparable to the central den- 
sity of atomic nuclei is made difficult by both the complexity of the interactions, 
discussed in the previous Chapter, and the approximations implied in any theoreti- 
cal description of quantum mechanical many-particle systems. 

The main problem associated with the use of the nuclear potential models de- 
scribed in Chapter [1] in a many-body calculation lies in the strong repulsive core of 
the NN force, which cannot be handled within standard perturbation theory. 

In non-relativistic many-body theory (NMBT), a nuclear system is seen as a 
collection of point-like protons and neutrons whose dynamics are described by the 
hamiltonian 



where t{i) and v{ij) denote the kinetic energy operator and the bare NN poten- 
tial, respectively, while the ellipses refer to the presence of additional many-body 
interactions (see Chapter [1]). 

Carrying out perturbation theory in the basis provided by the eigenstates of 
the noninteracting system requires a renormalization of the NN potential. This 
is the foundation of the widely employed approach developed by Briickner, Bethe 
and Goldstone, in which v{ij) is replaced by the well-behaved G-matrix, describing 
NN scattering in the nuclear medium (see, e.g. Ref . |40j ) . Alternatively, the many- 
body Schrodinger equation, with the hamiltonian of Eq. fl2.1l) . can be solved using 
either the variational method or stochastic techniques. These approaches have been 
successfully applied to the study of both light nuclei [25] and uniform neutron and 
nuclear matter [28ll26l[27lli2] . 

Our work has been carried out using a scheme, formally similar to standard per- 
turbation theory, in which nonperturbative effects due to the short-range repulsion 
are embodied in the basis functions. The details of this approach will be discussed 
in the following Sections. 




(2.1) 
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It has to be emphasized that within NMBT the interaction is completely deter- 
mined by the analysis of the exactly solvable two- and three-nucleon systems. As a 
consequence, the uncertainties associated with the dynamical model and the many- 
body calculations are decoupled, and the properties of nuclear systems ranging from 
deuteron to neutron stars can be obtained in a fully consistent fashion, without 
including any adjustable parameters. 



2.1 Correlated basis function theory 

The correlated states of nuclear matter are obtained from the Fermi gas (FG) states 
Inpc) through the transformation [43l l44] 

In the above equation, Inpc) is a determinant of single particle states describing 
noninteracting nucleons. The operator F, embodying the correlation structure 
induced by the NN interaction, is written in the form 

TV 

Fil,...,N)=S H (2.3) 

j>i=i 

where S is the symmetrization operator which takes care of the fact that, in general, 

[fijM ^ . (2.4) 

The structure of the two-body correlation functions fij must reflect the complexity 
of the NN potential. Hence, it is generally cast in the form (compare to Eq. fll.l7l) ) 

6 

^• = E/"M^^ (2-5) 

n=l 

with the O^- defined by Eq.f ll.l6|) . Note that the operators included in the above 
definition provide a fairly accurate description of the correlation structure of the 
two-nucleon bound state. The shape of the radial functions /"(^jj) is determined 
through functional minimization of the expectation value of the nuclear hamiltonian 
in the correlated ground state 

= {0\H\0) . (2.6) 

The correlated states defined in Eq. ( 12. 2p are not orthogonal to one another. However, 
they can be orthogonalized using an approach, based on standard techniques of 
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many-body theory, that preserves diagonal matrix elements of the hamiltonian 
Denoting the orthogonalized states by \n), the procedure of Ref. [15] amounts to 
defining a transformation T such that 

\n) \n) = f\n) , (2.7) 

with 

{n\H\n) = {n\H\n) . (2.8) 

Correlated basis function (CBF) perturbation theory is based on the decompo- 
sition of the nuclear hamiltonian 

H = Ho + Hj , (2.9) 

where Hq and Hj denote the diagonal and off-diagonal components of H, respec- 
tively, defined by the equations 

{m\HQ\n) = 6m.n{'m'\H\n) , (2.10) 
{m\Hi\n) = {l-6mn){m\H\n) . (2.11) 

The above definitions obviously imply that, if the correlated states have large over- 
laps with the eigenstates of H, the matrix elements of Hj are small, and the per- 
turbative expansions in powers of Hj is rapidly convergent. 

Let us consider, for example, the Green function describing the propagation of 
a nucleon in a hole state Il6] 



GiKuj) = (0|4 — — ^ ak|0)/(0|0) . (2.12) 

12 — rL(j — LU — 17] 

In the above equation, 77 = 0"'', aj^ and are creation and annihilation operators 
and the exact ground state |0), satisfying the Schrodinger equation H\0) = Eq\0) 
can be obtained from the expansion [ITJ HH] 



|0) = $^(-)"( 



Hi - AEq 
Ho - < 



|0) , (2.13) 



where AEq = Eq - E^ , with E^ defined by Eq.(l2J]). 

In principle, using Eq. fl2.13p and the similar expansion 

I = I yf.n( Hi-AEo y 

H-Eo-uj-ir] Ho- E^ -u -ir] ^ \Ho - E^ - u - ir] J ' ^ ' ^ 

the Green function can be consistently computed at any order in Hj. However, 
the calculation of the matrix elements of the hamiltonian appearing in Eqs. fl2.12p - 
fl2.14p involves prohibitive difficulties and requires the development of a suitable 
approximation scheme, to be discussed in the following Section. 
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2.2 Cluster expansion formalism 

The correlation operator of Eq.f l2.3l) is defined such that, if any subset of the parti- 
cles, say ii, . . .ip, is removed far from the remaining ip+i, . . . zat, it factorizes accord- 
ing to 

...,N)^ Fp{ii, . . . tp)FN-pitp+i, ...In) . (2.15) 

The above property is the basis of the cluster expansion formalism, that allows one 
to write the matrix element of a many-body operator between correlated states as 
a sum, whose terms correspond to contributions arising from isolated subsystems 
(clusters) involving an increasing number of particles. 

Note that in this Section we will use non normalized correlated states, defined 
as (compare to Eq. fl2.2p ) 

|n) = FIufg) . (2.16) 



2.2.1 Ground state expectation value of the hamiltonian 

Let us consider, as an example, the expectation value of the hamiltonian H in the 
correlated state |0), defined in Eq.f l2.16i) . We will closely follow the derivation of 
the corresponding cluster expansion given in Ref. [H] and neglect, for the sake of 
simplicity, the three body potential Vijk. Under this assumption, we can write the 
hamiltonian as in Eq. fl2.1l) . 

The starting point is the definition of the generalized normalization integral 

I{P) = {0\exp[PiH-n)]\0) , (2.17) 

where 

\p\<pf 

Pf being the Fermi momentum, is the ground state energy of the noninteracting 
Fermi gas at density p = 2p|,/37r^. Using the definition of Eq. fl2.17l) we can rewrite 
the expectation value of the hamiltonian in the form 



(OW) _ d 



(2.19) 

13=0 



The cluster property of F can be exploited to define a set of N\/{N — p)\p\ sub- 
normalization integrals, associated with each p-particle subsystem (p = 1, . . . ,A^) 

Ii{(3) = (^|exp[/3(t(l)-eO)]K) , 

/,,(/5) = (zj|4(12)exp[/3(t(l)+t(2) + t;(12)-e°-e°)]F2(12)|zj%, 
hMP) = m , (2-20) 
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where the indices i,j, . . . label states belonging to the Fermi sea, the ket l^i ■ ■ - in) 
describes n non interacting particles in the states ii . . . i„, e° = |p|f/2m is the kinetic 
energy eigenvalue associated with the state \i) and the subscript a refers to the fact 
that the corresponding state is antisymmetrized. For example, in the case of two 
particles 

\^j)a = ^ m - m) . (2.21) 

To express ln/(/3) in terms of the In Ii-^ i^{j3) , we start noting that 1^ is close to the 
product of li and Ij. It would be exactly equal if we could neglect the interaction, 
described by the potential f (12), and the correlations induced by both ^2(12) and 
Pauli exclusion principle. This observation suggests that lij can be written as 

lij I^IjYij , (2.22) 

with Yij ~ 1. Extending the same argument to the J's with more than two indices, 
we obtain 

h = 

I- = YYY- 

h...N = I = \{Y^ \{Y,, Y^ j, , (2.23) 



implying 



In /(/3) = J] In + ^ In y,,- + . . . + In Fi...^ . (2.24) 

i j>i 

It can be shown [H] that each term in the rhs of Eg. (12.241) goes like in the 
thermodynamic limit. In addition, the p-th term collects all contributions to the 
cluster development of ln/(/3) involving, in a connected manner, exactly p Fermi 
sea orbitals. Therefore, the p-th term can be referred to as the p-body cluster 
contribution to ln/(/5). 

The decomposition fl2.23p allows one to rewrite the expectation value of the 
hamiltonian in the form 



(H) = To + {AEh + (AE)3 + . . . + (A^);v (2.25) 

with 

. fi 

(2.26) 



dp 

ii<i2<...<ip 



/3=0 



Note that {AE)i = 0, as the above definitions imply 

h = Y, = l, (2.27) 
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and 

1^0, (2,28) 

To make the last step we have to use Eg. (12. 23 p and express {AE)p in terms of 
the Substitution of the resulting expressions 

= h 

Y,, = (2.29) 
: (2.30) 



into Eq. (l2:26|) with p = 2 yields 



1 dlii ( dli dij 



d(3 d(3 



(2.31) 

J/3=0 



The "normalizations" /ij|/3=o appearing in the denominator differ from unity by 
terms 0{1/N) at most, that can be disregarded in the N oo limit. As a result, 
we obtain 

(AE), = Wij , (2-32) 

where (see Eq.(l22Q])) 

= (U1 ^[^2(12), [t(l) +t(2), F2(12)] ] +F2(12)i;(12)F2(12) , (2.33) 

Note that in the above equation we have assumed that the correlation operator be 
hermitian, i.e. that ^2(12) = fJ(12) = /12 (see Eq.(l23D). The exphcit expression 
of (A£')2, in the case of six component potential and correlation operator, is given 
in Appendix [Bl 

Each term of the expansion fl2.25p can be represented by a diagram featuring p 
vertices, representing the nucleons in the cluster, connected by lines corresponding 
to dynamical and statistical correlations. The terms in the resulting diagrammatic 
expansion can be classified according to their topological structure, and selected 
classes of diagrams can be summed up to all orders solving a set of coupled integral 
equations, called Fermi hyper- netted chain (FHNC) equations [HI |50] . 



2.2.2 Transition matrix elements 

The nuclear matter response to an external probe delivering energy uj and momen- 
tum q can be written in the form 

^(q,a;) = |(/|O(q)|0)p5(c. + E, - Ef) , (2.34) 
/ 
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where O is the operator inducing a transition from the ground state |0), carrying 
energy Eq, to a final state |/), carrying energy Ej. In the simple case of interaction 
with a scalar probe, resulting in a density fiuctuation 

0(q) = p(q) = E«U«k, (2.35) 

k 

ajj. and Ok being nucleon creation and annihilation operators, respectively. 

In order to obtain the response within the CBF approach, the cluster expansion 
formalism discussed in the previous section must be extended to the case of transition 
matrix elements. 

Consider a (non normalized) correlated one particle-one hole final state 

|/) = |ph) = Fat,ah|O^G). (2.36) 
To obtain the response we need to calculate the matrix elements 



(Olp^lph) 



(0|0)l/2(ph|ph)V2 



(2.37) 



The cluster expansion of the above quantity can be carried out using a formalism 
somewhat different from the one described in the previous Section, originally devel- 
oped in Ref. [51]. The starting point is again the generalized normalization integral, 
that in this case is written in the form 



r ^ 

^ i=i 



■''^0„,(l)...0„^(iV), (2.38) 



In the above equation xi, . . . , xn denote the nucleon degrees of freedom, and 
$J(1, . . . ,A^) = (xi . . . xn\^fg) is the FG ground state wave function, i.e. the anti- 
symmetrized product of the single particle orbitals 0mi(l) • • • (pniN^i^)- Note that 
the calculation of matrix elements involving states describing Fermi systems, only re- 
quires the antisymmetrization of either the initial or the final state. The coordinate 
space expressions of the operators Pq and Wph are 

N N 
1=1 i=l 

i = l '^"^i^ ' 1 = 1 
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Acting on the product . . .0mjv5 replaces the hole state orbital with 
the particle state orbital (pp. In terms of generalized normalization integrals we can 
write 



{ph|pi|0) = /nd..*5FV'F*o=|^^5^ , (2-41) 



a=/3=0 



(ph|ph) = f l[dx,^l,^F^ F^ph = Wh(0,0) (2.42) 

where JP^'^'^(0,a) is obtained from Eq. (l2.38p replacing the FG ground state with the 
one particle-one hole state $p/i, and 

N 

(0|0) = / l[dxi^lF^ F% = /o,o(0,0) . (2.43) 

i=i 

To carry out the cluster expansion we need to define subnormalization integrals, 
involving an increasing number of orbitals 



= y t/xi0„,,(l)e^''^«e"^-'^«0„^(l) = X, 

= v^l rfxirfx2(/.L.^^,(l,2)4(12)e'^[^^(^)+''^(2)]^^(^2) 
X e°[^''^(^)+^^''(^)V™,(l)0™,(2) = X,Xj + Xi, 



h...N = l''''\P,a) = Xn-a ■ ■ ■ Xn...b , (2.44) 

where (^9^,^^.(1,2) = [0^, (1)0^^.(2) - 0.;„^(2)0.;„^,(l)]/y2), and the sum in the last 
line is extended to all partitions such that a + . . . + b = N. Note that in this 
case the (small) deviation between lij and the product lilj is characterized through 
their difference, rather than the ratio (compare to Eq. fl2.23l) ). 

The thermodynamic limit (i.e. the limit A,V oo, with A/V = const) of 
Ioph{P,o:) can be best identified rewriting it in the form originally obtained in Ref. 

m 

N 

Io,ph{P,a) = l[X,{P,a)e''o-''^^'''''^ , (2.45) 
1=1 

with 

Go,ph = J2^^^+OiN-') , (2.46) 

j>i 
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where 
Defining 

we finally obtain [51] 

(0|0)l/2(ph|ph)l/2 



' XiXj 



da 



-[Gph,pfe(0,0)+Go,o(0,0)]/2 



J a=0 



(3=0 



From the above equations it follows that, at two-body cluster level [51 



(Olp^lph) 



(0|0)l/2(ph|ph)V2 



E 



d dXi 



8(3 da 



f3=a=0 j>i 



dp da 



I3=a=0 



(2.47) 



(2.48) 



(2.49) 



(2.50) 



Note that the derivation of Eq. fl2.5Up has been carried out consistently with that 
of Eq. (12.311) . i.e. neglecting all contributions 0{N~^) to the normalization of the 
correlated two-nucleon states. 



2.3 Effective interaction 

At lowest order of CBF, the effective interaction V^s is defined by the equation 

(H) = {Ofg\To + V,s\Ofg) . (2.51) 

As the above equation suggests, the approach based on the effective interaction 
allows one to obtain any nuclear matter observables using perturbation theory in 
the FG basis. However, as discussed in the previous Section, the calculation of the 
hamiltonian expectation value in the correlated ground state, needed to extract V^g 
from Eq. (l2.5ip . involves severe difficulties. 

In this Thesis we follow the procedure developed in Refs. [52l[3l], whose authors 
derived the expectation value of the effective interaction by carrying out a cluster 
expansion of the rhs of Eq. (12.511) . and keeping only the two-body cluster contribu- 
tion. The resulting expression, that can be obtained from Eqs. (l2.26|l - (l2.33p through 
a simple rearrangement of the kinetic energy contributions, reads 



i<j 



i<j 

(2.52) 



--(VV12) - -(V/22) ■ V + v{12)fu 
m m 
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where the laplacian and the gradient operate on the relative coordinate. Note that 
feff defined by the above equation exhibits a momentum dependence due to the 
operator (V/jj) ■ V, yielding contributions to nuclear matter energy through the 
exchange termsQ However, our numerical calculations show that these contributions 
are small, compared to the ones associated with the momentum independent terms. 
As a consequence, the results presented in this Thesis have been obtained using only 
the static part of the effective interaction (12.521) . i.e. setting 

Vef,{^J) = k (-^ + v{^J)\ = v:^{n,)0- , (2.53) 

The properties of the operators O" with n = 1, . . . , 6, leading to the above result, 
are given in Appendix Rl 

The definition of fefr given by Eqs. (l2.52p and (12.531) obviously neglects the ef- 
fect of three-nucleon interactions, whose inclusion in the hamiltonian is known to 
be needed in order to explain the binding energies of the few-nucleon systems, as 
well as the saturation properties of nuclear matter. To circumvent this problem, 
we have used the approach originally proposed by Lagaris and Pandharipande [53] . 
in which the main effect of the three-body force is simulated through a density 
dependent modification of the two-nucleon potential at intermediate range, where 
two-pion exchange is believed to be the dominant interaction mechanism. Neglect- 
ing, for simplicity, the charge-symmetry breaking components of the interaction, the 
resulting potential can be written in the form 

v{^3)= [<(^«.)+<(r.,)e-^^'' + t;?(r.,)]Oi;. , (2.54) 

n=l,14 

where w", and denote the long- (one-pion-exchange), intermediate- and short- 
range part of the potential, respectively. The above modification results in a repul- 
sive contribution to the binding energy of nuclear matter. The authors of Ref. [53] 
also include the small additional attractive contribution 

^Etna = 72P'(3 - 2f3^)e-^'P , (2.55) 

with [3 = {Pp — Pn) / {Pp+ Pn) , where pp and p„ denote the proton and neutron density, 
respectively. The values of the parameters 71, 72 and 73 appearing in Eqs. (12.541) and 
(12.551) have been determined in such a way as to reproduce the binding energy and 
equilibrium density of nuclear matter [53] . 

Besides the bare two- and three body-potentials, the effective interaction is de- 
termined by the correlation operators fij defined by Eq. (l2.5p . The shapes of the 

^Thc direct contribution is vanishing, as it involves the integration of an odd function of k = 
ki — kj . 
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radial functions f^{Tij) are obtained from the functional minimization of the energy 
at the two-body cluster level, yielding a set of coupled differential equations to be 
solved with the boundary conditions 

^>--^^) = { J; : 12,3,4 ' (2.56) 
Fivij > rfj = , n = 5,6 (2.57) 

and 



dvij 



,n = 1,2,3,4 



drij 



, n = 5,6 , (2.58) 



rii=dt 



d and dt > d being variational parameters. The above conditions simply express 
the requirements that i) for relative distances larger than the interaction range the 
two-nucleon wave function reduces to the one describing non interacting particles 
and ii) tensor interactions have longer range. 

For any given value of nuclear matter density, we have solved the Euler-Lagrange 
equations resulting from the minimization of the binding energy at two-body cluster 
level, whose derivation and explicit form is given in Appendix [Cl using the values 
of d and dt reported in Ref . [51] . 

As an example, the results corresponding to nuclear matter at equilibrium density 
are illustrated in Fig J2.1l showing the central component of the correlation functions 
acting between a pair of nucleon carrying total spin and isospin S and T, respectively. 
The relations between the Jts of Fig J2.1l and the /" of Eq.f l2.5p are given in Appendix 
Rl The shapes of the J'ts clearly reflect the nature of the interaction. In the T = 
S = channel, in which the potential exhibits a strong repulsive core, the correlation 
function is very small at r < 0.5 fm. On the other hand, in the T = 5* = 1 
channel, the spin-isospin state corresponding to the deuteron, the repulsive core is 
much weaker and the potential becomes attractive at r > 0.7 fm. As a consequence, 
the correlation function does not approach zero as r — and exceeds unity at 
intermediate range. 

In Fig J2.2l the components of the effective interaction at equilibrium density are 
compared to the corresponding components of the truncated fg potential. It clearly 
appears that screening effects due to NN correlations lead to a significant quenching 
of the interaction. 

Figure [2^3] shows a comparison between the central [n = 1, left panel) and spin- 
isospin [n = 4, right panel) components of the effective interaction of Eq. fl2.53p . 
calculated at different densities using the Argonne fg potential. The density depen- 
dence is associated with the correlation functions, which depend on p through the 
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Figure 2.1. Interaction potentials (upper panel) and correlation functions (lower 
panel) acting in the spin-isospin channels 5 = and T = (solid lines) and S = 
and T = 1 (dashed lines). The potential is the Argonne and the correlation 
functions correspond to nuclear matter at equilibrium density. 



correlation ranges, d and dt, and the Fermi distributions. The jumps in the radial 
behavior of the effective interactions, clearly visible in Figs. 12.21 and 12.31 are due to 
the discontinuity in the second derivative of the correlation functions. 

2.3.1 Energy per particle of neutron and nuclear matter 

The effective interaction described in the previous Section was tested by computing 
the energy per particle of symmetric nuclear matter and pure neutron matter in first 
order perturbation theory using the FG basis. 
Let us consider nuclear matter at density 

4 

P = 5^ PA , (2.59) 

A=l 

where A = 1,2,3,4 labels spin-up protons, spin-down protons, spin-up neutrons and 
spin-down neutrons, respectively, the corresponding densities being px = x\p, with 
^^xa = 1. For example, for symmetric nuclear matter xi = 2:2 = = 0:4 = 1/4, 
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Figure 2.2. Comparison between the components of the bare Argonne poten- 
tial (dashed hnes) and the effective potential defined by Eq. (|2.53p (solid lines), 
calculated at nuclear matter equilibrium density. 
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Figure 2.3. Density dependence of the central (n = 1, left panel) and spin-isospin 
(n = 4, right panel) components of the effective interaction of Eq. (|2.53p . calculated 
using the Argonne fg potential. The dot-dash, dashed and solid lines correspond 
to p = 0.04, 0.16 and 0.32 fm~^, respectively. 



while for pure neutron matter Xi = X2 = and X3 = X4 = 1/2. Within our approach, 
the energy of such a system can be obtained from 

I = ^ E ^ + f E E ^A^^ / d'r v:, K - B:/{p,,,rnp^,,r)] . (2.60) 

A \^JL n 

In the above equation, ppx = (Gir'^pxY^^ and the Slater function i is defined as 

i{pF,xr) = Y^ e^''-'-^(pi.,A-|k|) . (2.61) 

k 

The explicit expression of the matrices 

A^^ = (A^|0"|A/i) , Bl = (A/i|0"|/iA) , (2.62) 

where |A/i) denotes the two-nucleon spin-isospin state, is given in Appendix lAl 

In Fig. 12.41 our results are compared to those of Refs. [H] and [27]. The calcula- 
tions of Ref. [H] (diamonds, results given in the sixth column of Table VI) have been 
carried out using a variational approach based on the FHNC-SOC formalism, with a 
hamiltonian including the Argonne f is NN potential and the Urbana IX three-body 
potential [38] • The results of Ref. [27] (dashed line of the lower panel) have been 
obtained using the Wg model and the same three-body potential, within the frame- 
work of the Auxiliary field diffusion Monte Carlo (AFDMC) technique. It appears 
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p [fm ] 



Figure 2.4. Energy per particle of symmetric nuclear matter (upper panel) and 
pure neutron matter (lower panel). The solid lines represent the results obtained 
using Eq. ()2.60p . whereas the diamonds correspond to the results of Akmal, Pand- 
haripande and Ravenhall [Jl]. The dashed line of the lower panel represents the 
results of the AFDMC approach or Ref . [27] . The square in the upper panel shows 
the empirical saturation point of symmetric nuclear matter. 



that the effective interaction approach provides a fairly reasonable description of the 
EOS over a broad density range. 

The comparison between the results of our calculation, based on the two-body 
cluster approximation, and those obtained taking into account higher order many- 
body effects deserves a comment. In view of the fact that the contribution of clusters 
involving more than two nucleons is known to be sizable, our approach has to be 
regarded as an effective theory, designed to provide lowest order results in agreement 
with the available "data". Effective theories are widely employed in many areas of 
Physics, including nuclear matter theory. For example, the Walecka model [55] is 
designed to reproduce the nuclear matter empirical saturation properties in the mean 
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field approximation, i.e. at tree level, although the corresponding loop corrections 
are known to be large. 

It is worth noting that the empirical equilibrium properties of symmetric nuclear 
matter are accounted for without including the somewhat ad hoc density dependent 
correction of Ref. [H]. The authors of Ref. [32] argued that this may be ascribed 
to the different description of the three-body force. It should also be emphasized 
that, using Vcs of Eg. (12.531) and the three-nucleon interaction (TNI) model of Ref. 
[53] , one effectively includes the contribution of clusters involving more than two 
nucleons. 

In addition to the correct binding energy per nucleon and equilibrium density 
{E/N = 15.96 MeV at p = 0.16 fm~^), our calculation also yields a quite reasonable 
value of the compressibility module, K = 230 MeV. 

It has to be kept in mind that our approach does not involve adjustable param- 
eters. The correlation ranges d and dt have been taken from Ref. [M], while the 
parameters entering the definition of the TNI have been determined by the authors 
of Ref. [53] through a fit of nuclear matter equilibrium properties. 



2.3.2 Single particle spectrum and effective mass 

Using the effective interaction, the effective mass can be obtained from the single- 
particle energies eA(p), that can be easily computed in Hartree-Fock approximation 
|46] . The resulting expression is (compare to Eq. (l2.60p ): 

""^^^ = 1^ + ^ E E / d'rv:^{r) [Al^ - i?,Xpr)^(p^,,r)] , (2.63) 

where p = |p| and jo is the spherical Bessel function: jo{x) = sin(x)/x. Figure 
12.51 shows e{p) = J2x^)^^>^{p) symmetric nuclear matter at equilibrium density. 
For comparison the corresponding results from Ref. [57] are also displayed. They 
have been obtained using the FHNC-SOC approach and the Urbana W14 two-nucleon 
potential, modified according to the TNI model of Lagaris and Pandharipande |53j . 
The nucleon effective mass, m*, is related to the single-particle energy through 

1 1 de / ,x 

— = - — . 2.64 

m* p dp 

The density dependence of the ratio m*{pF)/m of PNM, obtained from the Vcs 
discussed in this Chapter, is shown in Fig. 12.61 It is worth mentioning that for 
SNM at equilibrium, we find m^{j)p)/m = 0.65, in close agreement with the lowest 
order CBF result of Ref. [56] . The results of CBF calculations at second order show 
a ~20% enhancement of the effective mass at the Fermi surface, due to medium 
polarization effects [56] • We do not find this enhancement, as these effects are not 
taken into account in the Hartree-Fock approximation. 
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Figure 2.5. Solid Line: momentum dependence of the single particle energies 
obtained from the CBF effective interaction in the Hartree-Fock approximation 
(see Eq. ()2.63p ). The diamonds show the results of Ref . [57] . 
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Figure 2.6. Density dependence of the ratio m*'{pF)/'m of PNM, obtained from 
Eqs. (|2.63p and (|2.64l) using the effective interaction described in the text. 
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2.3.3 Spin susceptibility of neutron matter 

The results of numerical calculations show that the energy per particle of nuclear 
matter can be accurately approximated using the expression 



with 



a 
7 



ErP^ + E^rl^ 



(2.65) 



(Xg + X4) - (Xi +X2) 
(Xg - X4) - (Xi - X2) 



(2.66) 



In symmetric nuclear matter xa = 1/4 for all A (see the definition in Section r2.3.2p . 
yielding E/N = -Esnm = -^'o; while in pure neutron matter, corresponding to Xi = 
X2 = and X3 = X4 = 1/2, E/N = -EpNivi = Eq + Er, implying that Er can be 
identified with the symmetry energy. 

Let us consider fully spin-polarized neutron matter. The two degenerate states 
corresponding to X3 = 1 and X4 = (a = 1, spin- up) and X3 = and X4 = 1 
(a = — 1, spin-down) have energy. 



E^ = E^ = EpNM + E„ , (2.67) 
For arbitrary polarization a, the energy can be obtained from 



with Ecr = Ecr ■ 

the expansion 



E. 



E{a) = E{0) + 



dE 
da 



1 d^E 



a 



Q = 



2 da' 



a^ + 



(2.68) 



a=0 



As E must be an even function of a (see Eq. fl2.67p ). the linear term in the above 
series must be vanishing and, neglecting terms of order a^, we can write 



AE = E{a) - E{Q) = \ 



a 



(2.69) 



Q=0 



2 da^ 

In the presence of a uniform magnetic field B the energy of the system becomes 

Eeia) = E{a) - afiB, (2.70) 

where B denotes the magnitude of the external field, whose direction is chosen as 
spin quantization axis, and fi is the neutron magnetic moment. 
Assuming that equilibrium is achieved at a = ao, i.e. that 



dE 

da 



fiB = 



(2.71) 



a=aQ 
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Neutron Matter 




P/Po 



Figure 2.7. Ratio between the spin susceptibility obtained from Eqs. (j2.75p 
and (j2.60p and the FG model result. The points with error bars show the 
AFDMC results of Ref.ESI. 



we obtain 



«o = ( ^ ) • (2.72) 



a=0 



From the definitions of the total magnetization 



-1 



M = flips - p4) = l^aop = /^M ^ ) Bp, (2.73) 



a=0 



and the spin susceptibility x 
we finally obtain 



M = xB , (2.74) 



The above equation shows that, within our approach, the spin susceptibility of 
neutron matter can be easily calculated from Eq. fl2.60p 

Figure 12.71 shows the density dependence of the ratio between the susceptibility 
of neutron matter obtained from the effective interaction and that corresponding to 
the FG model. For comparison, the results of Ref . [58] , obtained within the AFDMC 
approach using the Argonne Vg NN potential and the Urbana IX three-body force. 
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are also displayed. It appears that the inclusion of interactions leads to a substantial 
decrease of the susceptibility over the whole density range, and that the agreement 
between the two theoretical calculations is remarkably good. 
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Chapter 3 

Nuclear matter response 



In this Chapter, we will discuss the response of nuclear matter to an external probe, 
defined in Eq. fl2.34p of Chapter [21 Extensive experimental studies of the nuclear 
response have been carried out mostly through inclusive electron scattering experi- 
ments (for a recent review see, e.g., [9]). The wealth of available data, corresponding 
to a variety of targets, ranging from the few nucleon systems, having A < 4, to nu- 
clei as heavy as Gold {A = 197) and Lead {A = 208), can be reliably extrapolated 
to the A oo limit to obtain quantitative empirical information on the nuclear 
matter response [59] . 

Electron scattering experiments have exposed the deficiencies of the independent 
particle model of nuclear dynamics. On the other hand, many body approaches 
explicitly including dynamical correlation effects provide a quantitative account of 
the measured cross sections in a broad kinematical domain [9j. 

As a pedagogical example, we will first consider the response to a scalar probe. 
The generalization to the case of electromagnetic and charged current weak inter- 
actions will be discussed in Chapters H] and [51 



3.1 Many-body theory of the nuclear response 

Within NMBT, the nuclear response to a scalar probe delivering momentum q and 
energy u, defined in Eq. (12.34]) . can be written in terms of the the imaginary part of 
the polarization propagator n(q,ci;) according to [IHl [H] 

S{q,u) = - Im n(q,a;) = - Im (0|pt ^ p^\0) , (3.1) 

TT 71 H — — UJ — 17] 

where ?7 = O"*" denotes an infinitesimal positive quantity and the operator pq, de- 
scribing the density fiuctuation induced by the probe, is given in Eq (l2.35p . 
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The above definition is best suited to estabhsh the relation between 5'(q,a;) and 
the nucleon Green function, leading to the popular expression of the response in 
terms of spectral functions [HI HTj. 

Equation (13.11) clearly shows that the interaction with the probe leads to a tran- 
sition of the struck nucleon from a hole state of momentum k, with |k| < pp, to 
a particle state of momentum k + q, with |k + q| > pp. Hence, the calculation of 
S'(q,a;) amounts to describing the propagation of a particle-hole pair through the 
nuclear medium. 

The Green function is the quantum mechanical amplitude associated with the 
propagation of a particle from x = (t,x) to x' = (t',x') [IH!- In uniform matter, due 
to translation invariance, it only depends on the difference x — x', and after Fourier 
transformation to the conjugate variable k = {k , E ) can be written in the form 

= Gh{k,E) + Gp{k,E) , (3.2) 

where Gh and Gp correspond to propagation of nucleons in hole and particle states, 
respectively. 

The connection between Green function and spectral functions is established 
through the Lehman representation [46] 

(3.3) 



Ph{k,E') Pp%E') 



E' -E-ir] E-E' -iri 



G{k,E) = j dE' 
implying 

Ph{k,E) = 5^ |(n(^,„i)(-k)|ak|0^)|25(E-E(-) + Eo) = ^ Im ^.(k,^) , (3.4) 

n 

Pp{k,E) = J]m^+i)(k)|ai|0,v)r5(E + EW -Eo) = ^Im Gp{p,E) , (3.5) 

n 

where |(?7,(Ar±i)(±k)) denotes an eigenstate of the {A ± l)-nucleon system, carrying 
momentum ±k and energy Eii^\ 

Within the FG model the matrix elements of the creation and annihilation op- 
erators reduce to step functions, and the Green function takes a very simple form. 
For example, for hole states we fined 



^Note that, according to our definitions, the hole spectral function is defined for E > — /i, /j, 
being the Fermi energy. 
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with = |k^|/2m, implying 

PfgA^^E) = 9{pp - \k\)6{E + el) . (3.7) 

Strong interactions modify the energy of a nucleon carrying momentum k according 
to el — y el + S(k,£'), where Il(k,E) is the comj)/ex nucleon self-energy, describing 
the effect of nuclear dynamics. As a consequence, the Green function for hole states 
becomes 

«'"'-^) - WTTlh^ ■ <^-**> 

A very convenient decomposition of G/i(k,£') can be obtained inserting a com- 
plete set of (A— l)-nucleon states (see Eqs. fl3.2p -f l3.4l) ) and isolating the contributions 
of one-hole bound states, whose weight is given by [60] 



Zk = K-k|ak|0)|2 = 9{pp - \k\)^k ■ (3.9) 



Note that in the FG model these are the only non-vanishing terms, and $fc = 1, 
while in the presence of interactions < 1- The resulting contribution to the Green 
function exhibits a pole at — e^, the quasi-particle energy being defined by the 
equation 

efe = e° + Re S(k,efc) . (3.10) 
The full Green function can be rewritten 

GhiKE) = ^.^\ . + Gf (k,E) , (3.11) 

E + Cfc + iZk Im S(k,efc) 

where is a smooth contribution, associated with (A — l)-nucleon states having 
at least one nucleon excited to the continuum (two hole-one particle, three hole-two 
particles . . . ) due to virtual scattering processes induced by nucleon-nucleon (NN) 
interactions. The corresponding spectral function is 

1 ZL Im S(k,ei,) ^r/-, „n 

vr + + Re S(k,efc)J^ + [Zfclm S(k,efc)J^ 

The first term in the right hand side of the above equation yields the spectrum 
of a system of independent quasi-particles, carrying momenta |k| < pp^ moving in 
a complex mean field whose real and imaginary parts determine the quasi-particle 
effective mass and lifetime, respectively. The presence of the second term is a con- 
sequence of nucleon-nucleon correlations, not taken into account in the mean field 
picture. Being the only one surviving at |k| > pp^ in the FG model this correlation 
term vanishes. 

Figure 13.11 illustrates the energy dependence of the hole spectral function of 
nuclear matter, calculated in Ref.[17| using CBF perturbation theory and a realistic 
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Figure 3.1. Energy dependence of the hole spectral function of nuclear matter. [17] 
The solid, dashed and dot-dash lines correspond to |k| = 1, 0.5 and 1.5 fm^^, respec- 
tively. The FG spectral function at |k| = 1 and 0.5 fm~^ is shown for comparison. 
The quasi-particle strengths of Eq. (j3.9p . are also reported. 



nuclear hamiltonian, including the Urbana Vi^ potential and the TNI discussed in 
the previous Chapter. Comparison with the FG model clearly shows that the effects 
of nuclear dynamics and NN correlations are large, resulting in a shift of the quasi- 
particle peaks, whose finite width becomes large for deeply-bound states with |k| <^ 
Pf- In addition, NN correlations are responsible for the appearance of strength at 
|k| > Pf- The energy integral 

n{k) = j dE Ph{k,E) (3.13) 

yields the occupation probability of the state of momentum k. The results of Fig. 
13.11 clearly show that in presence of correlations n{\k\ > pf) 7^ 0. 



3.2 Nuclear response and spectral functions 

In general, the calculation of the response requires the knowledge of Ph and Pp, as 
well as of the particle-hole effective interaction. |48[ [6T] The spectral functions are 
mostly affected by short range NN correlations (see Fig. 13.11) . while the inclusion 
of the effective interaction, e.g. within the framework of the Tamm-Dancoff ap- 
proximation (TD) or the Random Phase Approximation (RPA), [^I] is needed to 
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account for collective excitations induced by long range correlations, involving more 
than two nucleons. 

At large momentum transfer, as the space resolution of the probe becomes small 
compared to the average NN separation distance, ^(q,^;) is no longer significantly 
affected by long range correlations. The authors of Ref. [62] found that for |q| > 500 
MeV RPA corrections are negligibly small, if computed using finite size interactions. 

In this kinematical regime the zero-th order approximation in the effective inter- 
action, according to which hole and particle propagate independent of one another, 
is expected to be applicable. The corresponding response can be written in the 
simple form 



The widely employed plane wave impulse approximation (lA) [9] can be readily 
obtained from the above definition replacing Pp with the FG result, which amounts 
to disregarding final state interactions (FSI) between the struck nucleon and the 
spectator particles. The resulting expression reads 



^.^(q,^) = / d'kdE Phik,E)9{\k + q| - Pf)S{u; -E- ej^+q,) . (3.15) 



Figure [3^ showing the uj dependence of the nuclear matter structure function at 
|q| = 5 fm~^, illustrates the role of correlations in the target ground state. The solid 
and dashed lines have been obtained from Eg. 03.151) using the spectral function of 
Ref. [17] and that resulting from the FG model (shifted in such a way as to account 
for nuclear matter binding energy), respectively. It clearly appears that the inclusion 
of correlations produces a significant shift of the strength towards larger values of 
energy transfer. 

At moderate momentum transfer, both the full response and the particle and 
hole spectral functions can be obtained using non relativistic many-body theory. 
The results of Ref. [17] suggest that the zero-th order approximations of Eqs. fl3.14l) 
and (13.151) are fairly accurate at |q| > 500 MeV. However, in this kinematical 
regime the motion of the struck nucleon in the final state can no longer be described 
using the non relativistic formalism. While at lA level this problem can be easily 
circumvented replacing the non relativistic kinetic energy with its relativistic coun- 
terpart, including the effects of FSI in the response of Eq. (l3.14p involves further 
approximations, needed to obtain the particle spectral function at large |q|. 

3.3 Particle spectral function at large momentum 

A systematic scheme to include corrections to Eq. (l3.15p and take into account FSI, 
originally proposed in Ref. [HI] , is discussed in Ref. [nS| • The main effects of FSI on 




(3.14) 
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Figure 3.2. Nuclear matter SiAi^,^) (see Eq. (j3.15p ). as a function of ui at |q| = 5 
fm~^. The solid and dashed lines correspond to the spectral function of Ref. [47j 
and to the FG model (shifted in such a way as to account for nuclear matter binding 
energy), respectively (taken from Ref. [63j ) . 

the response are i) a shift in energy, due to the mean field of the spectator nucleons 
and ii) a redistributions of the strength, due to the coupling of the one particle-one 
hole final state to n particle-n hole final states. 

In the simplest implementation of the approach of Refs.[64l EH], the response is 
obtained from the lA result according to 

^(q,^;) = J du' SiA{q.,uj')fq{u - u') , (3.16) 

the folding function /q being related to the particle spectral function through 

Pp(k + q,c^ -E) = e{kF - |k + q|) /|k+q|(c^ -E- efk+q,) (3.17) 



with e|k+q| = v^|k + q^ + m^. In the absence of FSI, /q shrinks to a (5-function and 
the lA result of Eq. (13.151) is recovered. 

Obviously, at large q the calculation of Pp(k + q,uj — E) cannot be carried out 
using a nuclear potential model. However, it can be obtained form the measured 
NN scattering amplitude within the eikonal approximation. The resulting folding 
function is the Fourier transform of the Green function describing the propagation 
of the struck particle, travelling in the direction of the z-axis with constant velocity 
v: 

f Juj) = [ — e-^e^/o'^i'M.iK) . (3.18) 
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where k + q q and 



Mq|W = (0|^Er|q|(r..+z)|0) 

j>i 



(3.19) 



In the above equation, r|q| is the Fourier transform of the NN scattering amph- 
tude at incident momentum |q| and momentum transfer |t|, A|q|(fc), parameterized 
according to 



(3.20) 



In principle, the total cross section a, the slope (3 and the ratio between the real 
and the imaginary part, a, can be extracted from NN scattering data. However, the 
modifications of the scattering amplitude due to the presence of the nuclear medium 
are known to be sizable, and must be taken into account. The calculation of these 
corrections within the framework of NMBT is discussed in Ref . |66j . 

In Eg. (13. 191) . the expectation value is evaluated in the correlated ground state. 
It turns out that NN correlation, whose effect on Ph is illustrated in Fig. 13.11 also 
affect the particle spectral function and, as a consequence, the folding function of 
Eq. (13.171) . Neglecting all correlations 



|q| 



|q| 



-vpa{i + a) 



(3.21) 



and the quasi-particle approximation 



Pp(q,u; - E) 



UJ 



E 



Im 



|q| 



+ 



Im 



(3.22) 



is recovered. 

Correlations induce strong density fluctuations, preventing two nucleon from 
coming close to one another. The joint probability of finding two particles at posi- 
tions ri and r2 can be written 



p(ri,r2) = {0\Y,6iv, - ri)5(r, - r2)|0) = p'gi\r^ - r^\) 

j>i 



(3.23) 



The above equation defines the radial distribution function g{r), which describes 
correlation effects. Figure 13.31 shows the typical shape of the radial distribution 
function resulting from the CBF calculation of Ref. [67] . 

The effect of correlation on FSI can be easily understood keeping in mind that 
the response is only sensitive to rescattering taking place within a distance ~ l/|q| 
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Figure 3.3. Radial distribution function of nuclear matter at equilibrium density, 
obtained from CBF perturbation theory using a realistic hamiltonian |67j . 

of the primary interaction vertex As the probabihty of finding a spectator within 
the range of the repulsive core of the NN force (< 1 fm) is small, the probability that 
the struck particle rescatter against one of the spectators within a length ~ l/|q| 
is also very small at large |q|. Hence, inclusion of correlations leads to a significant 
suppression of FSI. 

Fig. 13.41 shows the u dependence of the nuclear matter response of Eq. (13.161) 
at |q| =5 fm~^. The solid and dashed lines have been obtained using the spectral 
function of Ref.[17], with and without inclusion of FSI according to the formalism 
of Ref . [64] . respectively. For reference, the results of the FG model are also shown 
by the dot-dash line. The two effects of FSI, energy shift and redistribution of the 
strength from the region of the peak to the tails, clearly show up in the comparison 
between solid and dashed lines. 




^Note that this is no longer true in the case in which the hadronic final state is also observed. 
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Figure 3.4. Nuclear matter 5(q,u;), defined in Eq. (l3.16p . as a function of lj at 
|q| = 5 fm^^. The solid and dashed lines have been obtained from the spectral 
function of Ref.[17|, with and without inclusion of FSI, respectively. The dot-dash 
line corresponds to the FG model (shifted in such a way as to account for nuclear 
matter binding energy), respectively (taken from Ref. [63]). 
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Chapter 4 

Impulse Approximation regime 



As pointed out in the previous Chapter, the nuclear response has been extensively 
investigated carrying out inclusive electron scattering experiments. 

The first attempts to provide a quantitative estimate the of the measured electron- 
nucleus cross section were based on oversimplified models of nuclear dynamics. At 
the end of the seventies, Moniz suggested that the target may be described as a 
degenerate gas of protons and neutrons at given constant density p [68], the effect 
of the interactions being crudely taken into account by an average binding energy 
e. Despite its simplicity, the FG model of Ref. [M] was able to give a fairly accu- 
rate account of the electron-nucleus cross section in the region of the quasi-elastic 
peak, corresponding io xb = Q'^/2muj ~ 1, where is the Bjorken variable, 

= IqP ~ ^"^1 q, uj and m denote the momentum and energy transfer and the 
nucleon mass, respectively. 

In the past twenty years, with the availability of new data, extending in the region 
of high |q| and low cj, corresponding to xb ^ 1, the limits of the FG model, and 
more generally of all independent particle models, became apparent. Away from 
the quasi elastic peak correlation effects, not included in the mean field picture, 
become more and more important and the FG model is not longer able to describe 
the measured cross section. 

The experimental investigation of the neutrino-nucleus cross section involves 
additional difficulties due to the low counting rates and the lack of neutrino beams of 
fully specified properties. However, a quantitative understanding of the weak nuclear 
response is needed in a variety of different fields, ranging from nuclear astrophysics 
to the analysis of neutrino oscillation experiments. 

Electron scattering data provide a stringent test for validation of theoretical 
models of the nuclear response, in view of their application to the case of weakly 
interacting probes. For example, the success of the FG model in explaining electron 
scattering in the quasi elastic region at |q| < 500 MeV prompted its extension to 
neutrino scattering 
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In this Chapter we will review the application of the formalism on NMBT to the 
calculation of the electromagnetic and charged current weak cross sections in the 
region of large momentum transfer, where the lA is expected to be safely applicable. 



4.1 Electron- nucleus cross section 

The differential cross section of the process 

e + A^e' + X, (4.1) 

in which an electron carrying initial four- momentum k = (-E^k) scatters off a nu- 
clear target to a state of four-momentum k' = [E'^,\i'), the target final state being 
undetected, can be written in Born approximation as (see, e.g., Ref. [70] ) 

(fa _ a" E', j^u.yA .42^ 

where a is the fine structure constant. The leptonic tensor, that can be written, 
neglecting the lepton mass, as 

L^''' = 2[k,k', + k,k'^-g^,{kk')\ , (4.3) 

is completely determined by electron kinematics, whereas the nuclear tensor W^^, 
contains all the information on target structure. Its definition involves the initial and 
final hadronic states |0) and carrying four-momenta andpx, respectively, as 
well as the nuclear electromagnetic current operator J'^: 

W^^ = Y,{^\J^\X){X\rm^'\p, + q-px) , (4.4) 

X 

where the sum includes all hadronic final states. Comparison with Eg. ( 12.341) shows 
that the above tensor is the generalization of the nuclear response to the case of 
vector interaction. 

Calculations of at moderate momentum transfers (|q| < 0.5 GeV) can be 
carried out within nuclear many-body theory (NMBT), using non-relativistic wave 
functions to describe the initial and final states and expanding the current operator 
in powers of |q|/m, m being the nucleon mass (see, e.g., Ref. [71], [721 [73]). On the 
other hand, at higher values of |q|, corresponding to beam energies larger than ~ 1 
GeV, the description of the final states \X) in terms of non-relativistic nucleons is 
no longer accurate. Calculations of in this regime require a set of simplifying 
assumptions, allowing one to take into account the relativistic motion of final state 
particles carrying momenta ~ q as well as the occurrence of inelastic processes, 
leading to the appearance of hadrons other than protons and neutrons. 
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4.2 The impulse approximation 

As stated in Chapter [3|, the main assumptions underlying the impulse approximation 
(lA) scheme are that i) as the spatial resolution of a probe delivering momentum q is 
~ l/|q|, at large enough |q| the target nucleus is seen by the probe as a collection of 
individual nucleons and ii) the particles produced at the interaction vertex and the 
recoiling (A— l)-nucleon system evolve independently of one another, which amounts 
to neglecting both statistical correlations due to Pauli blocking and dynamical Final 
State Interactions (FSI), i.e. rescattering processes driven by strong interactions. 

In the lA regime the scattering process off a nuclear target reduces to the in- 
coherent sum of elementary processes involving only one nucleon, as schematically 
illustrated in Fig. 14. 1[ 




2 



Figure 4.1. Schematic representation of the lA scheme, in which the nuclear 
cross section is replaced by the incoherent sum of cross sections describing 
scattering off individual bound nucleons, the recoiling (A — l)-nucleon system 
acting as a spectator. 



Within this picture, the nuclear current can be written as a sum of one-body 
currents 

J'^T.J"^ (4-5) 

i 

while the final state reduces to the direct product of the hadronic state produced at 
the electromagnetic vertex, carrying momentum p^, and the (A — 1) -nucleon residual 
system, carrying momentum pji = q — Px (for simplicity, we omit spin indices) 

\X) ^ |x,p,) ® |7^,p7^) . (4.6) 

Using Eq. (14.61) we can rewrite the sum in Eq. (14.41) replacing 

\X){X\ I d^Px\x,Vx){VxA 5^rfV|7^,P7^)(P7^,7^| . (4.7) 

Substitution of Eqs. (I4.5I) -( 14T1) into Eq. (14. 4p and insertion of a complete set of free 
nucleon states, satisfying 

J £p\N,p){p,N\ = I , (4.8) 
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results in the factorization of the current matrix element 

1/2 



(0| J^X) = ^ (0|7^,p7^; N, - p^^) J](-p7e,iV|jf |x,p,.) , (4.9) 

leading to 

^mI' = 5^ /c^'P7^rf'px|(0|7^,p7^;N,-p^'■^ 



X 



J] (-p^,N|jf |a;,P.)(p.,x|jT|N, - Ptz) (4.10) 
X 5(=^)(q-p7^-p^.)<^(^ + ^o-^7^-^x), 



where Ep^ = \/\pTz^~-\-m? . Finally, using the identity 

S{uj + Eq-Eti-E^) = j dE5{E - m + Eq - E-n) 6{u - E + m - E^) , (4.11) 

and the definition of the target spectral function given in the previous Chaptei0, 

P(p,E) = 5^ |(0|7^, - p; N,p) |2 S{E-m + Eo- En) , (4.12) 



n 

we can rewrite Eq.f l4.4p in the form 



= Y.Jd'p (^) ^(P'^) ' (4.13) 



with Ep = Y Ip^I + m? and 

•^^-5^(p,N|jf|a;,p + q)(p + q,a;|jT|N,p) 5{uj + + - E,) . (4.14) 



Note that the factor (m/ a/Ipt^P + m?Y/'^ in Eg. (14. 9 p takes into account the implicit 
covariant normalization of (— p7^,A^| in the matrix element of jf. 

The quantity defined in the above equation is the tensor describing electromag- 
netic interactions of the i-th nucleon in free space. Hence, Eq. f l4.14p shows that in 
the lA scheme the effect of nuclear binding of the struck nucleon is accounted for 
by the replacement 

g = (^,q) ^ g = (2;,q) , (4.15) 



"'^As we will consider a target having N = Z = A/2, the spectral functions describing proton 
and neutron removal will be assumed to be the same. 
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with (see Eqs. and gUD) 



UJ 



UJ 
UJ 



Eq — E-ji 
E + m — 



(4.16) 



in the argument of u'f'^. This procedure essentially amounts to assuming that: i) a 
fraction 6uj of the energy transfer goes into excitation energy of the spectator system 
and ii) the elementary scattering process can be described as if it took place in free 
space with energy transfer uj = u — 6uj. This interpretation emerges most naturally 
in the |p| <C m limit, in which Eq. (14.161) yields Suj = E. 

Collecting together all the above results we can finally rewrite the doubly differ- 
ential nuclear cross section in the form 



da 



lA 



dVtpidEpi 



d^p dE P{p,E) 



^ daep 



+ {A-Z) 



dap. 



dVLpidEf.! 



dde'dEgi 
6{uj — E + m — Ex), 



(4.17) 



where daeN/dfle'dEe' {N = n,p denotes a proton or a neutron) is the cross section 
describing the elementary scattering process 



e{k) + N{p) e'{k') +x{p + q) 



given by 



da, 



eN 



dQe'dEp' 

stripped of both the flux factor and the energy conserving (5-function. 



(4.18) 



(4.19) 



4.3 The nuclear spectral function 

Non-relativistic NMBT provides a fully consistent computational framework that 
has been employed to obtain the spectral functions of the few-nucleon systems, 
having A= 3 [TH [75], [76] and 4 [TTJ [IS], as well as of nuclear matter, i.e. in the 
limit A ^ oo with Z=A/2 [17], [HD]. Calculations based on G-matrix perturbation 
theory have also been carried out for oxygen [8T1 [82] . 

The spectral functions of different nuclei, ranging from Carbon to Gold, have 
been modeled using the Local Density Approximation (LDA) |83], in which the ex- 
perimental information obtained from nucleon knock-out measurements is combined 
with the results of theoretical calculations of the nuclear matter P[p,E) at different 
densities [55] . 
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Nucleon removal from shell model states has been extensively studied by coinci- 
dence (e,e'p) experiments (see, e.g., Ref. [22]). The corresponding measured spectral 
function is usually parameterized in the factorized form 

Pmf(p,^) = \Mp)\^KiE - En) , (4.20) 

n 

where (pnip) is the momentum-space wave function of the single particle shell mode 
state n (e.g. Woods-Saxon wave functions), whose energy width is described by the 
function Fn{E — En) (e.g. a lorentzian). The normalization of the n-th state is given 
by the so called spectroscopic factor Z„ < 1, and the sum in Eq. fl4.20p is extended 
to all occupied states. Typically, Pmf{p,E) vanishes at E larger than ~ 30 MeV 
and IpI larger than ~ 250 MeV. Note that in absence of NN correlations the full 
spectral function could be written as in Eq. (14.201) . with Fn{E — En) = S{E — En) 
and Zn = 1. 

Strong dynamical NN correlations give rise to virtual scattering processes leading 
to the excitation of the participating nucleons to states of energy larger than the 
Fermi energy, thus depleting the shell model states within the Fermi sea. As a 
consequence, the spectral function associated with nucleons belonging to correlated 
pairs extends to the region of \p\ ^ pp and E ^ ep, where ep denotes the Fermi 
energy, typically < 30 MeV. 

The correlation contribution to P{p,E) of uniform nuclear matter has been cal- 
culated by Benhar et al for a wide range of density values [83]. Within the LDA 
scheme, the results of Ref. [83] can be used to obtain the corresponding quantity 
for a finite nucleus of mass number A from 

P^oAp.E) = J d'r pA{r)PZ!^, (p,E; p = p^(r)) , (4.21) 

where Pa(i") is the nuclear density distribution and P^^ {p,E; p) is the correlation 
component of the spectral function of uniform nuclear matter at density p. 
Finally, the full LDA nuclear spectral function can be written 

PLDAip,E) = Pmf{P,E) + P,orr{p,E) , (4.22) 

the spectroscopic factors Zn of Eq. fl4.20p being constrained by the normalization 
requirement 

J d^p dE Plda{p.E) = 1 . (4.23) 

The LDA spectral function of ^^O obtained combining the nuclear matter results 
of Ref. [83] and the Saclay {e,e'p) data [S^ is shown in Fig. 14. 2[ The shell model 
contribution Pmf{p,E) accounts for ~ 80 % of its normalization, whereas the re- 
maining ~ 20 % of the strength, accounted for by Pcorr{p,E), is located at high 
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Figure 4.2. Three-dimensional plot (left panel) and scatter plot (right 
panel) of the oxygen spectral function obtained using the LDA approxima- 
tion described in the text. 

momentum (|p| ^ pp) and large removal energy {E ^ ep). It has to be empha- 
sized that large E and large p are strongly correlated. For example, ~ 50 % of the 
strength at |p| = 320 MeV is located at E > 80 MeV. 

The LDA scheme rests on the premise that short range nuclear dynamics is 
unaffected by surface and shell effects. The validity of this assumption is confirmed 
by theoretical calculations of the nucleon momentum distribution, defined as 

n(p) = jdE P{p,E) = (0|at,ap|0) , (4.24) 

where and Cp denote the creation and annihilation operators of a nucleon of 
momentum p. The results clearly show that for A> 4 the quantity n{p)/A be- 
comes nearly independent of A in the region of large |p| (> 300 MeV), where NN 
correlations dominate (see, e.g., Ref. [85]). 

In Fig. I4.3l the nucleon momentum distribution of ^^O, obtained from Eq. (14.241) 
using the LDA spectral function of Fig. 14.21 is compared to the one resulting from 
a Monte Carlo calculation |86], carried out using the definition of Eq. (I4.24p and a 
highly realistic many-body wave function |88]. For reference, the FG model momen- 
tum distribution corresponding to Fermi momentum pp = 221 MeV, currently used 
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in the analysis of neutrino oscillation experiments (see, e.g. Ref.[87]), is also shown 
by the dashed line. It clearly appears that the n(p) obtained from the spectral 
function is close to that of Ref . [86] , while the FG distribution exhibits a completely 
different behavior. 




12 3 4 



P [fm-^] 

Figure 4.3. Momentum distribution of nucleons in the oxygen ground state. 
Solid line: LDA approximation. Dashed line: FG model with Fermi momentum 
Pf = 221 MeV. Diamonds: Monte Carlo calculation carried out by S.C. Pieper 
[H6] using the wave function of Ref. [88] . 

A direct measurement of the correlation component of the spectral function of 
^^C, obtained measuring the {e,e'p) cross section at missing momentum and energy 
up to ~ 800 MeV and ~ 200 MeV, respectively, has been recently carried out 
at Jefferson Lab by the E97-006 Collaboration [89]. The data resulting from the 
preliminary analysis appear to be consistent with the theoretical predictions based 
on LDA. 

4.4 Comparison to electron scattering data 

We have employed the formalism described in the previous Sections to compute the 
inclusive electron scattering cross section off oxygen at 0.2 <Q^< 0.6 GeV^ [901I9T]. 

The lA cross section has been obtained using the LDA spectral function shown 
in Fig. 14.21 and the nucleon tensor defined by Eq. fl4.14p . that can be written as 

< - < {-r . f:) . < (p^ - f (p" - f,^) . (4.25, 
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where p = (i?p,k) and the off-shell four momentum transfer q is defined by Eqs. 
(14.151) and (I4.16p . The two structure functions Wi and W2 are extracted from 
electron-proton and electron-deuteron scattering data. In the case of quasi-elastic 
scattering they are simply related to the electric and magnetic nucleon form factors, 
and Gm^^ through 



Numerical calculations have been carried out using the Hohler-Brash parameteri- 
zation of the form factors [92l [93], resulting from a fit which includes the recent 
Jefferson Lab data |94j . 

In the kinematical region under discussion, inelastic processes, mainly quasi-free 
A resonance production, are also known to play a role. To include these contribu- 
tions in the calculation of the inclusive cross section, we have adopted the Bodek 
and Ritchie parameterization of the proton and neutron structure functions |95j . 
covering both the resonance and deep inelastic region. 

In Figs. I4.4M.7I the results of our calculations are compared to the data of Ref. 
[97] . corresponding to beam energies 700, 880, 1080 and 1200 MeV and electron 
scattering angle 32°. For reference, the results of the FG model corresponding to 
Fermi momentum pp = 225 MeV and average removal energy e = 25 MeV are also 
shown. The results including FSI effects have been obtained from the approach 
described in Chapter [3l using the gaussian parameterization of Eq. (l3.20p . with the 
parameter values resulting from the fit of Ref. [96] 

Overall, the approach described in the previous Sections, involving no adjustable 
parameters, provides a fairly accurate account of the measured cross sections in the 
region of the quasi-free peak. On the other hand, the FG model, while yielding a 
reasonable description at beam energies 1080 and 1200 MeV, largely overestimates 
the data at lower energies. The discrepancy at the top of the quasi-elastic peak 
turns out to be ~ 25 % and ~ 50 % at 880 and 700 MeV, respectively. 

The results of NMBT and FG model also turn out to be sizably different in the 
dip region, on the right hand side of the quasi-elastic peak, while the discrepancies 
become less pronounced at the A-production peak. However, it clearly appears that, 
independent of the employed approach and beam energy, theoretical results signif- 
icantly underestimate the data at energy transfer larger than the pion production 
threshold. 

In view of the fact that the quasi-elastic peak is correctly reproduced (within an 
accuracy of ~ 10 %), the failure of NMBT to reproduce the data at larger uj may be 
ascribed to deficiencies in the description of the elementary electron-nucleon cross 




(4.26) 




(4.27) 
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Figure 4.4. Cross section of the process ^^0{e,e') at beam energy 700 MeV and 
electron scattering angle 32°. Solid line: full calculation, with inclusion of final 
state interaction effects. Dot-dash line: lA calculation, carried out neglecting 
FSI effects. Dashed line: FG model with pp = 225 MeV and e = 25 MeV. The 
experimental data are from Ref . [97| . 




Figure 4.5. Same as in Fig. 14. 4^ but for beam energy 880 MeV. 



section. In fact, as illustrated in Fig. 14. 8[ the calculation of the lA cross section at 
the quasi-elastic and A production peak involves integrations of P{p,E) extending 
over regions of the {p,E) plane almost exactly overlapping one another. 

To gauge the uncertainty associated with the description of the nucleon structure 
functions and w^, we have compared the electron-proton cross sections obtained 
from the model of Ref. [HS] to the ones obtained from the model developed in Refs. 
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Figure 4.6. Same as in Fig. 14.41 but for beam energy 1080 MeV. 
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Figure 4.7. Same as in Fig. 14.41 but for beam energy 1200 MeV. 

[98| [99| 1100] and from a global fit [lOlj including recent Jefferson Lab data [102j . 
The results of Fig. KB show that at = 1200 MeV and 9 = 32° the discrepancy 
between the different models is not large, being ~ 15 % at the A production peak. 
It has to be noticed, however, that the models of Refs. [95l EHl [9H1 [IDQl [M] have 
all been obtained fitting data taken at electron beam energies larger than 2 GeV, so 
that their use in the kinematical regime discussed in this work involves a degree of 
extrapolation. 

On the other hand, the results obtained using the approach described in this 
paper and the nucleon structure functions of Ref. [HS] are in excellent agreement 
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Figure 4.8. The solid and dashed hnes enclose the integration regions in the 
{p,E) plane relevant to the calculation of the lA cross section at the top of 
the quasi-elastic and A production peak, respectively, for beam energy 1200 
MeV and scattering angle 32°. 




a [GeV] 

Figure 4.9. Cross section of the process e + N ^ e' + X above pion production 
threshold, at beam energy 1200 MeV and scattering angle 32°. Solid line: H2 fit of 
Ref. [lOOj for ep scattering; dashed line: fit of Ref. [95] for ep scattering; diamonds: 
fit of Ref. |101j for ep scattering; dot-dash line: fit of Ref. [95] for en scattering. 

with the measured (e,e') cross sections at beam energies of few GeV [83] . 

Figure 14.91 also shows the prediction of the Bodek and Ritchie fit for the neutron 
cross section, which turns out to be much smaller than the proton one. The results of 
Ref. |103j suggest that extrapolating the Bodek and Ritchie fit to the low region 
relevant to tha analysis of the data of Ref. may lead to sizably underestimate the 
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Figure 4.10. lA cross section of the process ^^0(e,e') at beam energy 1200 
MeV and scattering angle 32°. Dashed hne: quasi-elastic; dots: quasi- 
free A production; dashes: nonresonant background; solid line: total. The 
experimental data are from Ref. |97j . 



neutron contributions. On the other hand, the fit of Ref. [95] consistently includes 
both resonant and nonresonant contributions to the nuclear cross section. In this 
regard, it has to be pointed out that the nonresonant background is not negligible. 
As illustrated in Fig. 14.101 for beam energy 1200 MeV and scattering angle 32° it 
provides ~ 25 % of the cross section at energy transfer corresponding to the A peak. 



4.5 Neutrino-nucleus cross section 

The Born approximation cross section of the weak charged current process 

lye + A^r +X , (4.28) 
can be written in the form (compare to Eq. (14.21) ) 

where G = Gp cos 6c, Gp and 9c being Fermi's coupling constant and Cabibbo's 
angle, Ee is the energy of the final state lepton and k and k' are the neutrino and 
charged lepton momenta, respectively. Compared to the corresponding quantities 
appearing in Eq. (14. 2p . the tensors L^^, and W'^'^ include additional terms resulting 
from the presence of axial- vector components in the leptonic and hadronic currents 
(see, e.g., Ref. [TUij l 
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Within the lA scheme, the cross section of Eq. (I4.29P can be cast in a form similar 
to that obtained for the case of electron-nucleus scattering (see Eq. (14.171) ). Hence, 
its calculation requires the nuclear spectral function and the tensor describing the 
weak charged current interaction of a free nucleon, w'^. In the case of quasi-elastic 
scattering, neglecting the contribution associated with the pseudoscalar form factor 
Fp, the latter can be written in terms of the nucleon Dirac and Pauli form factors 
Fi and F2, related to the measured electric and magnetic form factors Ge and Gm 
through 




(4.30) 
(4.31) 
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Figure 4.11. Differential cross section da/dQedi^ for neutrino energy = 1 GeV 
and electron scattering angle 9e = 30°. The lA results are represented by the 
dashed line, while the solid line corresponds to the full calculation, including the 
effects of FSI. The dotted line shows the prediction of the FG model with Fermi 
momentum kp = 225 MeV and average separation energy e = 25 MeV. 

Figure HTTT] shows the calculated cross section of the process ^^0{h'e,e), corre- 
sponding to neutrino energy = 1 GeV and electron scattering angle 6^ = 30°, 
plotted as a function of the energy transfer p = E^ — Ef>. Numerical results have been 
obtained using the spectral function of Fig. 14.21 and the dipole parameterization for 
the form factors, with an axial mass of 1.03 GeV. 

Comparison between the solid and dashed lines shows that the inclusion of FSI 
results in a sizable redistribution of the lA strength, leading to a quenching of the 
quasi-elastic peak and to the enhancement of the tails. For reference, we also show 
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the cross section predicted by the FG model with Fermi momentum Pf = 225 MeV 
and average separation energy e = 25 MeV. Nuclear dynamics, neglected in the 
oversimplified picture in terms of noninteracting nucleons, clearly appears to play a 
relevant role. 




[GeV^] 



Figure 4.12. Differential cross section da/dQ'^ for neutrino energy E = 1 GeV. 
The dot-dash line shows the I A results, while the solid and dashed lines have been 
obtained using the modified spectral function of Eq. (j4.32p . with and without 
inclusion of FSI, respectively. . 



It has to be pointed out that the approach described in Chapter [3|, while includ- 
ing dynamical correlations in the final state, does not take into account statistical 
correlations, leading to Pauli blocking of the phase space available to the knocked- 
out nucleon. 

A rather crude prescription to estimate the effect of Pauli blocking amounts to 
modifying the spectral function through the replacement 

P{p,E) Pip,E)ei\p + q| - pp) (4.32) 

where pp is the average nuclear Fermi momentum, defined as 

Pf = j p^(rK(r), (4.33) 

with Pf{'^) = (3vr^pyi(r)/2)^/^, Pa{y) being the nuclear density distribution. For 
oxygen, Eq. (14.331) yields pp = 209 MeV. Note that, unlike the spectral function, 
the quantity defined in Eq. (14.32^ does not describe intrinsic properties of the target 
only, as it depends explicitly on the momentum transfer. 



61 



4 - Impulse Approximation regime 



> 



o 



W 

T3 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 
[GeV] 



> 



o 



zi. 

w 

T3 



3 
2.5 

2 
1.5 

1 

0.5 




16, 



+ O ^ |J. + X 
Ey = 2 GeV 




0.5 



1.5 



E^ [GeV] 



Figure 4.13. Quasi-elastic differential cross section da/dE^ as a function of the 
scattered energy for the neutrino energy E = 0.8 and 2.0 GeV. The sohd line 
shows lA calculation with Pauli blocking as in Eq. (40), the dot-dash line lA 
calculation without Pauli blocking, and the dashed line FG model. 



The effect of Pauli blocking is hardly visible in the differential cross section shown 
in Fig. I4.1H as the kinematical setup corresponds to > 0.2 GeV^ at the quasi- 
elastic peak. The same is true for the electron scattering cross sections discussed in 
the previous Section. On the other hand, this effect becomes very large at lower Q^. 

Figure H:.12l shows the calculated differential cross section da/dQ"^ for neutrino 
energy = 1 GeV. The dashed and dot-dash lines correspond to the lA results 
with and without inclusion of Pauli blocking, respectively. It clearly appears that 
the effect of Fermi statistic in suppressing scattering shows up at < 0.2 GeV^ 
and becomes very large at lower Q^. The results of the full calculation, in which 
dynamical FSI are also included, are displayed as a full line. The results of Fig. 14.121 
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suggest that Pauli blocking and FSI may explain the deficit of the measured cross 
section at low with respect to the predictions of Monte Carlo simulations |105j . 

Figure 14.131 shows the i^^-nucleus cross sections as a function of the scattered 
muon energy, by comparing the cross sections calculated by FG, and by the use of 
the spectral function with and without Pauli blocking. Figure 14.131 shows that FG 
yields a larger high-energy peak contribution than the other two. This is not due to 
the Pauli blocking, but due to the nuclear correlation effects in the spectral function: 
the muons tend to be scattered with a higher energy. This effect should show up in 
the forward angle cross section and may have a direct effect on neutrino oscillation 
measurements. 
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Chapter 5 

Low momentum transfer regime 



In this chapter we focus on the nuclear matter response to weak interactions in the 
regime of low momentum transfer (|q| ~ 10 MeV), where the non relativistic ap- 
proximation is expected to be applicable. Within this approach, the initial and final 
states can be obtained from NMBT, while the weak current entering the definition 
of the tensor W^'^ (see Eg. (14.41) and (14.291) ) is expanded in powers of |q|/?Ti. At 
leading order, the resulting response can be written in the simple form (compare to 
Eq.dlMD) 

Sici,uj) = ^J2(^\Ol\n){n\0^\0)6iuj + Eo-En) . (5.1) 

n 

where, in the case of charged current interactions, Oq is the operator corresponding 
to Fermi or Gamow- Teller transitions. 

The calculation of the nuclear matter response has been carried out using CBF 
states, obtained from the states of the noninteracting system trough the transforma- 
tion (12. 2p . and the two-body cluster approximation for the weak transition matrix 
elements (see Chapter [2]). The effect of long range correlations, which are known to 
play a critical role at low momentum transfer, has been also included in our scheme 
in a fully consistent fashion, using the effective interaction defined in Chapter [2] and 
the Tamm-Dancoff (TD) approximation 1106] . We will restrict our discussion to 
the case of Fermi transitions. The extension to Gamow- Teller transitions is trivial. 



5.1 Non relativistic reduction of the weak charged 
current 

The starting point of our calculation is the non relativistic reduction of the weak 
charged current operator. Basically, one needs to expand in powers of |q|/m the 
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matrix elements 

(pVmIp) = u{p\s')T^AT^u{p,s) r]l,T+r]t , (5.2) 

where 

r'' = YF, + la^'q^.p- + q^Fs , (5.3) 

with a^^ = z[7^,7j,]/2, 7^ being the Dirac matrices, and 

= ^''^^Fa + q^l^'Fp + z7V'^"g,FT . (5.4) 

In the above equations, 7^ = Z7''7^7^7^, u(p,s) is the spinor describing a Dirac 
fermion with momentum p and spin polarization s, rjt is the Pauli spinor specifying 
the isospin state of the nucleon, r"*" = (r^ + ir'^)/'! is the isospin raising operator 
and the form factors Fi, F2, F^, Fa-, Fp and Ft are functions of the squared four 
momentum transfer q^. 
From the definition 

u{p,s) = N^i J , (5.5) 

where A^p is a normalization constant and Xs is a Pauli spinor, it follows that, to 
zero-th order in |q|/m, we can write 

(pVoy(O)lp) ^ i^i(O) x\.Xs rilr+T], . (5.6) 

Due to the antidiagonal structure of the matrix 7^, the corresponding matrix el- 
ement of the axial vector current, (p'| 7(^(0) |p), does not have any zero-th order 
contributions. 

In conclusion, denoting gy = Fi{0) the non relativistic reduction amounts to 
making the replacement 

(r^l Jj-V^(q)|r,) ^ Of (q) = gvS{r, - e^^^rt , (5.7) 

which defines the Fermi transition operator Of (q) in coordinate space. 

The axial part of the current contributes through the yU. = i = 1,2,3 components 

(p'|Ja|p) = m(p',s')[Fa775 + ^p7oq75]M(p,s 77t/r+?7t 

= m^(p',s')[^a7o775 + ^p7oq75]r^M(p,s) 77t'r+?7t . (5.8) 

In the above equation, the zero-th order term proportional to Fp vanishes because 
7o75 is antidiagonal. As for the term proportional to Fa we find instead 

W7.^(i \)( ° ^V? iVf^ !V (5.) 



-/ y V -cr y V / y v o o- 
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Finally, the vector part of J+ does not contribute at zero-th order. 
Making use of Eq. (15.91) , we can then write 

(r:|J+(q)|r,) ^ Or(q) = gA^ir, - e^^'^^a^r^ , (5.10) 

with qa = Fa{0), which defines the operator inducing Gamow- Teller transitions. 

5.2 Correlated matrix elements and effective op- 
erators 

Using correlated states implies severe difficulties in the explicit calculation of the 
weak matrix element. In the FG model, the nuclear response is non vanishing 
only when the final nuclear state differs from the initial state for the presence of a 
particle excited outside the Fermi sea and a hole in the Fermi sea. In the presence of 
correlations, which can induce virtual nucleon-nucleon scattering processes leading 
to excitation of nucleons to states outside the Fermi sea, more complex scenarios 
must also be considered. For example, if the initial state has a two particle-two hole 
component, the final state can be a three particle-three hole state or, if the probe 
interacts with an excited nucleon, a two particle-two hole state. 

In the following we will consider only the dominant transition, between the cor- 
related ground state and a correlated one particle-one hole (ph) state. The corre- 
sponding Weak matrix element can be written 

M,. = <'"'l^'f^l°> . . (5.11) 

where F is the correlation operator defined in Eq. (l2.3p . Here the kets |0) and \ph) 
correspond to the ground and one particle-one hole Fermi Gas states, respectively, 
and O = ^ - Oi, Oi being the Fermi or the Gamow- Teller transition operator (see 
Eqs.dEZD) and flSlOD l 

In calculating the weak matrix element, we will use the two-body cluster approx- 
imation discussed in Chapter O Let us define 

9ij = fij - 1 , (5.12) 

with fij defined as in Eq. (12.51) . Note the the gij is short ranged, and therefore its 
matrix elements are small. 

At two-body level, the cluster expansion Eq. (l5.11l) yields 



{ph\F^OF\0) ^ {ph\{l + Y,9^,)0{l + Y.9.m 



j>i j>i 
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N(N - 1) 

N{ph\Oi\0) + ^-{ph\{Oi + O2,9i2}\0) 



N{N -1) 



{ph\gi2{O, + O2)gi2\0) 



(5.13) 



where {A,B} = AB + BA. The above equation, suggest the definition of an effective 
operator Oll^ , acting on Fermi Gas states, reminiscent of the effective interaction of 
Chapter [2l From 

{ph\F^OF\0) = {ph\0'f^O). 
it follows that, at the two-body cluster level (compare to Eq.f l5.13i) ) 



(5.14) 



1 

N' 



iV- 1 



{Oi + 02.912} 



N 



■[guiOi + 02)912] 



(5.15) 



Note that the O^^^ is a two-body operator, as it includes screening effects arising 
from nucleon-nucleon correlations. 

Replacing Eq. fl5.13l) into Eq. fl5.1ip leads to 



Mph -N1 + N2 + N3 



where 



N2 



N{ph\Oi\0) 
N{N-1) 



N{N-1) 



{ph\{Oi + O2,gi2}\0) , 
{ph\gi2iOi + O2)gi2\0) 



(5.16) 

(5.17) 
(5.18) 

(5.19) 



The wave function of the FG ground state \I^o(-R) = (-R|0) can be written as a 
determinant according to 



^o(^) 



1 



0i(l) 



(Pn{N) 



(5.20) 



where (pmi'n) is the wave function describing the n-th particle in the state m, with 
momentum and spin and isospin projections Sm and tm, respectively 



One particle-one hole states, and the corresponding wave functions '^ph{R) 
can be built from the ground state through 

\ph) = aW\Q) 



(5.21) 
(Ob/i), 

(5.22) 
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using Eq. lKWh . 

As Qij is a two-body operator, it turns out to be convenient rewriting the ground 
and one particle-one hole wave functions, \l/o(-R) and in the form (see 

Eq.(I52DD) 

MR) = ^^(^ _ X^(-l)""+"^0.(l)0/3(2)1>,^.,,3(3, ■ ■ ■ ,N), (5.23) 

%h{R) = _ I 5^ Yl (-1)""^"^ 0„(1)0^(2)$,,,^„,^,,(3, ■ ■ ■ ,iV) 

+ 5^(-l)"'^+"'^0„(l)0,(2)$,^„,,(3---iV) 

+ ^(-l)""+"'^+Va(2)0p(l)$,^„,,(3, ■ ■ ■ iV) 

where $p^-y^Q^/3(3, ■ ■ ■ ,A^) denotes the wave function of a — 2-particle system, with 
a particle in the state p and holes in the state 'j ^ a, (3. 

The (N-2)-particle wave functions satisfy the following orthonormalization rela- 
tions 

(5.24) 



j d'r, ■ ■ ■ $i^«,/3(3, ■ ■ ■ ,A^)$p,-,^,.,;3',h(3, ■ ■ ■ ,iV) = 



rfVg ■ ■ ■ rf\„<l>:^^„^(3, ■ ■ ■ ,iV)$^^„.,;,(3, ■ ■ ■ ,iV) = 6^^, V . (5.25) 

With the help of the above equations, the weak matrix element can be reduced to 
an integral over the coordinates of two particles. 

The next section will be devoted to the details of the calculation of the Fermi 
transition matrix element. We will start writing it as a sum of contributions cor- 
responding to the physical processes schematically represented by the diagrams of 
Figs. I5.m5.5[ The different contributions will be labeled by three indices 

• the first index, B or C, denotes the order (first and second, respectively) in 
the correlation g; 

• the value of the second index, 1 or 2, indicates that the Fermi operator acts on 
particle 1 (any nucleon inside the Fermi-sea in the initial state) or on particle 
2 (the active particle, carrying momentum h and spin projection ah in the 
ground state); 

• the third index specifies the direct (d) and exchange (e) contributions to the 
matrix elements. Note that exchange terms carry an additional minus sign. 
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As the Fermi operator involves the isospin raising operator, the only non van- 
ishing matrix elements are those in which the particle excited outside the Fermi-sea 
is a proton and the active particle in the initial state is a neutron. 

5.2.1 Fermi Transition 

We will now discuss the explicit calculation of the weak matrix element. Since the 
corresponding operator does not induce spin transitions, all contributions to the 
matrix element involve a Kronecker delta expressing the condition ap = at, where 
(Tp and ah are the spin projections of the particle and hole states, respectively. In 
order to simplify the notation, this factor will be omitted. 
From Eq. fl5.17p we obtain 

= p (27r)=^(5(=^)(p-h-q) , (5.26) 
where Xp{h) = Xa,,,^ and r]p^h) = Vr.^.y 




X 

Figure 5.1. Diagram associated with the contribution A'^i (see Eg. (15. 171) ). 

The above equation shows that A^^i, the leading term of the expansion, is just 
the Fermi transition matrix element in the absence of correlations. 

lo order to write down the complete expression of 7V2, let us consider the matrix 
element 

Bl = ^^i^Lzll j dml^{0,,g,,]^>,{R) (5.27) 
= f c?'nf/V2[0„(l)0p(2)-0p(l)0„(2)]{Oi,(?i2}0„(l)0,.(2) , 

a 

in which we have used Eqs.f l5.23l) -( l5.24p and the sum is extended over the Fermi 
sea. The right hand side of Eq. fl5.28p has both direct and exchange contributions; 
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for the direct term we can write 
1 1 

X {ap\{Tt,On\ah) 



Bid 



a n 



(5.28) 



where the Qni'r') are the radial functions entering the definition of and r = |ri— r2|. 
From now on, the ket indicates a state in which particles 1 and 2 have spin- 
isospin a and /3, respectively. Carrying out the change of variables 



(5.29) 



r = n - r2 

[ R=(ri + r2)/2, 
and integrating over R we obtain 

Slrf = £p(2vr)35(3)(p_q_h) d're^''^gn{r){ap\{Tt,Ol,'}\ah), (5.30) 



N 



E 



(5.31) 



where we have used 

E^ E 

Let us now focus on the spin-isospin matrix element. From the definition of the 
six operators {\1.16\\ and 



{a\T~^\a) 







''^^{ap\cri(T2\ah) = , 

'^^{ap\Si2\ah) = , 

act 

it follows that the only non- vanishing contribution is given by 

X] {ap\{T^,TiT2}\ah) = 8 . 



(5.32) 



(5.33) 



Note that the previous result was obtained exploiting the anticommutator 

{ri+,TiT2} = 2t^. (5.34) 
The final expression of Bid turns out to be 

Bid = ^{2n)^S^^\p - q - h) 8 j d^re'"^' g2{r). (5.35) 
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Figure 5.2. Diagram associated with the contribution Bid (see Eq. (l5.28p ). 



Let us now consider the exchange term. Performing again the change of variables 
of Eq. fl5.29p . we find the expression 



3re~^i^+^">gr.ir){pa\{T^,0l'}\ah) , (5.36) 



i?le = ii^(2vrf5(3)(p-q-h)5^ / d 

a,n 

which, after summing over ka, becomes 

Ble = ^(27r)35(3)(p - q - h) J d're~'''^ gnir)iikpr){pa\{T+ ,0'^}\ah) 



where 



pJo (2vr)3 (kprr 
is the Slater function. 

For the spin-isospin matrix element, using the completeness relation 



(5.37) 
(5.38) 



we can write: 



J](pa|{l,r+}|a/i) = 2 
^ (pa|{riT2,ri+}|a/i) = 2 
^ (pa|{cricr2,r+}|a/i) = 6 
{pa\{TiT2cricr2,T^}\ah) = 6 



(5.39) 
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{pa\{Si2,T^}\ah) = , 
^ {pa\{Si2TiT2,r^}\ah) = . 



Collecting all the above results we finally obtain 



Ble 



_P_ 

AV 



(27r)35(3) (p _ q - h) j d^re-'^^iikpr) [2^i + 2^2 + 6^3 + 6^4] • (5.40) 



Note that, in order to simplify the notation, in the above equation and in the rest 



Figure 5.3. Diagram associated with the contribution Ble (see Eq. (l5.40p ). 

of this Section the functional dependence of the QnS on r is omitted. 

The contributions B2d and B2e can be obtained from Bid and Ble through the 
substitution 



For the direct term, B2d, we find 



B2d = ^(27r)35(3)(p -q-h) 8 j d^r gi( 



(5.41) 



(5.42) 



while the exchange term, B2e, reads 

P 




Figure 5.4. Diagram associated with the contribution B2d (see Eq. (j5.42p l. 



B2e = ^(27r)35(3)(p -q-h)J d\e-''^H{kFr)[2g^ + 2^72 + 6(73 + 6^74] . (5.43) 
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The contributions of Eq. fl5.19p . N3, can again be obtained from the corresponding 
contributions of Eg. (15. 171) . Ni, through the substitution 

Y,9n{a'P'\{0l',rt}\aP) ^ Y.9n9n'{a'P'\OnrtO'J\aP) . (5.44) 



As the six operators O^^ form an algebra (see Appendix , we can write 

0l.'0l^ = J2^r.mr0l' . (5.45) 



Exploiting this property, the product of two or more operators can be rewritten as 
a linear combination of the six operators. We will also make use of the following 
relations 



(^A)(^B) 

(cricr2)^ 
Si2{cricr2) 

r.2 
^12 



(AB) + zo-(A X B), 
3 - 2(cricr2), 
{cricr2)Si2 = Si2, 
6 + 2(crio-2) - 25i2. 



(5.46) 
(5.47) 
(5.48) 
(5.49) 



The first equation of the above group follows from the properties of the Pauli matri- 
ces. The second can be obtained directly from the first, while the third and fourth 
can be easily derived writing the tensor S12 in the form: 



= ^(Sfifj - 5ij)alal, 



(5.50) 



where f is the unit vector in the direction of r. Using Eqs. fl5.46p -f l5.49p we find 




Figure 5.5. Diagram associated with the contribution Cld (see Eq. fl5.51]) ). 



Cld = ^(27r)35(3)(p_q_h) 

X J d're'^'[8gig2 + 8gl + 24gl + A8gl + 2Ag^g^ + A8g,g,] , (5.51) 
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and 

Cle = ^(2nf6^'\p-q-h) J d're-'''^£ikprM + 2g,g2 + Qgm + Q9m 

+ gl - 'igl - 'igl + + 12(?6 + %2'?3 + Qg2gi - Qg^g^ + lAg^g^] . (5.52) 

Following the same procedure and using the substitution (I5.4ip we also obtain 



C2d = ^(27r)35(3)(p_q_h) j d'r[8gl-4gl + 12gl-12gl + 24gl-24g 



(5.53) 



C2e = ^(27r)35(3)(p_q_h)y" d're-''^'[2gl + 2g,g2 + 6g^g, + 6gig, 

+ gl - + 12^5' + 12^2 + Qg^g,_^ _^ Q^^^^ _ Qg^g^ ^ 24^5^0] " (5.54) 

5.3 Calculation of the response 

In the previous Sections, we have defined all the ingredients needed to obtain the 
weak response of nuclear matter in the one particle-one hole channel. For reasons 
that will become apparent in the next Section, devoted to the discussion of the TD 
approximation, numerical calculations have been carried out in a cubic box of finite 
volume L^, using a discrete set of one particle one-hole states. Within this scheme, 
the allowed nucleon momenta k are given by 

+ + (5.55) 

Ij 

where n^, n^, = 0, ±1, ±2, ... and x, y and z are unit vectors along the 
directions of the cartesian axes. 

As we are interested in the low momentum transfer regime (|q| <^ kp)-, we 
determine the size of the normalization box by requiring that q be on the lattice, 
i.e. that 

27r 

q = — (rigX + n^y + n^^) . (5.56) 

For fixed momentum transfer, i.e. for fixed n^, n^, and |q|, the above equation 
yields 

L = 27rj^ , (5.57) 

with 



'^. = VK)' + K)' + K)'- (5-58) 

Obviously, in the L — s> oo limit our procedure must reproduce the results obtained 
working with a continuum set of one particle-one hole states. 
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Let us rewrite the definition of tlie response using tlie effective operator, defined 
as in Eq. fl5.14p . associated with Fermi transitions 

Sp{q,u;) = ^ J2 \{ph\0!ff{q)\0)\'6{uj + Eo-E,,)6{p-h-c[) . (5.59) 

h|<fc^,|p|>fef' 

Note that the factor 2 in the right hand side of the above equation comes from the 
sum over the particle and hole spin projections. 

For the sake of illustration, we will first consider the noninteracting FG model. 
In this case we find 

|M,,(q)r = \{ph\O!f^{q)\0)\' - \{ph\0^ {com' = 1 , (5-60) 

and 

Eph - Eo = ujp - ujh = ^ph , (5.61) 

with 

IkP 

a;^ = 1^ . (5.62) 
2m 

The sum appearing in Eq. fl5.59p is extended to all lattice momenta 

27r 

hi = — (nfx + rviy + <£) (5.63) 

such that |hj| < kp and |hj + q| > kp, with q given by Eq. fl5.56p . 

Obviously, the response obtained from the discrete set of final states consists of 
a collection of delta function peaks located at o; = u;j = ci^ihi+qi — ojhi- A smooth 
function of u has been obtained using a finite width gaussian representation of the 
energy conserving 5 function, i.e. replacing 

5{u - u{) exp{-[(cu - u.;)/af} . (5.64) 

(Ta/TT 

For sufficiently small values of cr, the results obtained using this procedure become 
independent of a. 

In Fig. 15.61 the FG response at |q| = 0.3 fm~^, obtained from an analytical 
calculation using a continuum set of one particle-one hole states, is compared to 
that resulting from the above procedure with (n^, n^, n^) = (1,2,3). This choice 
corresponds to a normalization box of linear dimension L = 78 fm, containing ~ 
78400 nucleons. The corresponding number of one particle-one hole states in the 
basis is ~ 3000. It clearly appears that the basis is large enough to reproduce the 
results obtained in the continuum limit. 
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0.010 




5 10 15 20 

0) [MeV] 

Figure 5.6. Comparison between the FG response at |q| = 0.3 fm^^ evaluated 
in a cubic box of side L = 78 fm using a basis of ~ 3000 states (diamonds) 
and that resulting from an analytical calculation using a continuum set of one 
particle-one hole states (solid line). 

5.4 Correlated Hatree-Fock approximation 

The inclusion of interaction leads to sizable modifications of the FG response. Cor- 
relation effects in the transition matrix elements, taken into account through the 
use of the effective operator, produce a quenching of the Fermi transition matrix 
elements of ~ 15 %, largely independent of the hole momentum, as illustrated in 
Fig. 15.71 Note that the effect on the response is larger, as its definition involves the 
transition probabilities, i.e. the squared matrix elements. This feature is apparent 
in Fig. 15.81 where the FG response, represented by diamonds, is compared to that 
obtained using the effective operator in the calculation of Mph and the FG spectrum 
of Eq. (15.621) . represented by crosses. 

An even larger modification is produced by interaction effects on the single parti- 
cle energies. Corrections to the kinetic energy spectrum (I5.62p have been calculated 
within the Hartree-Fock (HF) approximation, using the CBF effective interaction 
discussed in Chapter [2] (see Eq.( ]2.63l) ). Replacing the single particle energies of 
Eq. (l5.62i) with the HF energies, shown by the solid line in Fig. 12.51 leads to a 
sizable broadening of u region corresponding to non-vanishing response. 

The squares of Fig. 15.81 show the response evaluated within the correlated HF 
approximation , i.e. using the CBF effective operator and the HF spectrum ob- 
tained from the CBF effective interaction, at |q| = 0.3 fm"'^. Comparison with the 
diamonds and the crosses shows that using the HF spectrum leads to an increase of 
the upper limit of the energy transfer u from ~ 18 MeV to ~ 27 MeV. 
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Figure 5.7. Fermi transition matrix element at |q| = 0.3 fm~^^ as a function of 
the magnitude of hole momentum |h|. The dashed horizontal line corresponds 
to the result of the FG model. 
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Figure 5.8. Nuclear matter response at |q| = 0.3 fm~ . Diamonds: FG 
model. Crosses: results obtained using the effective operator in the calcu- 
lation of Mph and the FG spectrum of Eq. (15.62p . Squares: correlated HF 
approximation described in the text. 



In Figs. I5.9tl5.10l the comparison between FG model and correlated HF approxi- 
mation is extended to larger values of the momentum transfer |q| = 1.8, and 3.0 fm~^, 
corresponding to the regions |q| > kp and |q| > 2kF, respectively. In the case 
|q| = 1.8, fm~^ we have used a box of side L =~ 47 fm, corresponding to a basis 
of ~ 3900 states, while the calculation at |q| = 3.0, fm~^ has been carried out with 
L =~ 44 fm and ~ 3300 basis states. From Figs. I5.8H5.10I it clearly appears that 
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the quenching of the peak and the shift of the strength towards larger values of uj, 
resulting from the inclusion of interaction effects, are sizable in all instances. 
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Figure 5.9. Nuclear matter response at |q| = 1.8 fm~^. Diamonds and 
squares correspond to the FG model and the correlated HF approximation 
described in the text. 
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Figure 5.10. Same as in Fig. 15.91 for momentum transfer |q| = 3.0 fm ^. 



5.5 Tamm-Dancoff Approximation 

In the previous section we have discussed the nuclear response in the correlated 
HF approximation, in which the bare Fermi transition operator is replaced by the 
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effective operator of Eq. fl5.14p and the final state is assumed to be a one particle-one 
hole state. 

It is important to realize that the FG one particle-one hole states, while being 
eigenstates of the HF hamiltonian, defined as 

Hhf = X^ej , (5.65) 

i 

with Cj given by Eg. (12. 631) . are not eigenstates of the full nuclear hamiltonian. As a 
consequence, there is a residual interaction Vres that can induce transitions between 
different one particle-one hole states, as long as their total momentum q, spin and 
isospin are conserved. 

In order to include the effects of these transitions, we use the TD approximation, 
which amounts to expanding the final state in the basis of one particle-one hole states 
according to jl06j 

I/) = |q, TSM) = J2cr''\K P^ = h, + q, TSM) , (5.66) 

i 

where S and T denote the total spin and isospin of the particle hole pair and M is 
the spin projection. 

At fixed q, the excitation energy of the state (15.661) . u^, is obtained solving the 
eigenvalue equation 

H\f) = iEo + iuf)\f) , (5.67) 

where Eq is the ground state energy. Substituting Eq. (15.661) into Eq. (l5.67p and 
multiplying by {hj, pj, TSM\ from the left leads to 

J2{h„P„TSM\H\h,,p,,TSM)cr'' = 0^"" 

= (Eo + ^^^^^cf ^ , (5.68) 

with 

Hjl^^' = {Eo + ep, - eh,)S,, + (h„ p„ T5M|Ke.|h„ p„ TSMs) , (5.69) 



and p] 

(h,-, p,-, TSM\Vres\K p., TSM) = {v,ff)o - {Veff)E • (5.70) 

In the above equation, Vejf is the CFB effective interaction discussed in Chapter [2l 
and the direct and exchange contributions to the matrix elements are given by 

{veff)D = ^J dh,d?r^ e-'ip.r^+i^^r,) ^a.^a^M'ef f\aH,ap^) e'^^^^-^+P'^^) , (5.71) 
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(''eff)E =y^j d\,dh, e-^(h'^^+P^^^) {ar,^a,XfMh,^P.) e^^^^^-^+P^^^^ , (5.72) 

where the two- nucleoli state \ah(y-p) describes a hole in the spin-isospin state ah 
and a particle in the spin-isospin state a^, coupled in such a way as to obtain the 
assigned values of T and S and M. 

We can now exploit the fact that, in the case of Fermi transitions, the requirement 
that the matrix elements of Of^j be non- vanishing implies that the particle hole pair 
be in a 5* = 1 state with M = ±1. As a consequence, H^j^^ , cf^'^ and uf^'^^ become 
independent of M. 

For any fixed T and S, the diagonalization of the hamiltonian matrix defined in 
Eqs. fl5.69p - (l5.72p determines the eigenvalues u"^^ and the the corresponding eigen- 
vectors c^^ with n = 1, 2, . . . , Nb, Nb being the number of states in the one 
particle-one hole basis. 

Collecting all the above results, we can finally write the response in the TD 
approximation as 



2 

T=0,1 M=±l n 



5^(c^i).(h„ p„ TlM|Oe//(q)|0) 



2 

Tl 



5{uj-ui') , (5.73) 



where (c^^), denotes the i-th. component of the eigenvector belonging to the eigen- 
value u'^^. 

The main features of the response in the TD approximation can be understood 
considering a simple model in which Vejf is assumed to be central and spin-isospin 
independent, and the exchange contribution to Eq. (15.701) is neglected. As a result 
we can make the replacement 

(h„ p„ T5M|K.e.|h„ p,, TSM) -^^^ j dh veff{r) e'^' = VefM) , (5.74) 

leading to the eigenvalue equation 

The above equation can be solved graphically plotting the right hand side as a 
function of uj and finding the intersections with the line J^{uj) = 1, as shown in Fig. 
15.111 for |q| = 0.3 fm~^ and a basis consisting of 8 states. The upper and lower 
panel correspond to Veff{q) = —.06 and 0.7 MeV, respectively. It appears that in 
the latter case the spectrum exhibits an eigenvalue lying well outside the particle 
hole continuum, corresponding to a collective excitation, reminiscent to the plasmon 
mode of the electron gas. 
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Figure 5.11. Graphical solution of Eg. (15. 751) for |q| = 0.3 fm~^ and a basis 
consisting of 8 states. The upper and lower panel correspond to Veff{q) = —0.06 
and —0.7 MeV, respectively. 



The nuclear matter response at |q| = 0.3 fm~^ obtained using the TD approx- 
imation is shown in Fig. I5.12[ The solid line in the lower left panel corresponds 
to the full calculation, including both the direct and exchange contributions to the 
matrix elements of the effective interactions (see Eq. fl5.70p ). while the solid line in 
the upper left panel has been obtained neglecting the exchange part. In both left 
panels, the dashed line refers to the correlated HF approximation. All calculation 
have been carried out using ~ 3000 basis states. 

The full TD response exhibits a sharp isolated peak corresponding to the collec- 
tive mode, lying ~ 4 MeV above the upper limit of the particle hole continuum. If 
only direct contributions are included, the peak is still clearly visible, but located 
at lower u, and not well separated from the continuum. 

The right panels of Fig. 15.121 illustrate the contributions to the TD response 
associated with different final states, corresponding to total isospin of the particle 
hole pair T = (solid lines) and T = 1 (dashed lines). In the case of Fermi 
transitions, the total spin is = 1 and the two allowed spin projections, M = ±1 
give equal contributions. The right lower and upper panels show the results of the 
full calculations and those obtained neglecting the exchange terms. 
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Figure 5.12. Nuclear matter response in the TD approximation at |q| = 0.3 fm~^. 
The soHd Hne in the upper left panel corresponds to the full calculation, while that 
in the lower left panel has been obtained neglecting the exchange term in Eq. (j5.70p . 
The dashed lines show the results of the correlated HF approximation. The right 
panels illustrate the contributions of the T = (solid lines) and T = 1 (dash-dot 
lines) to the TD response, evaluated including direct and exchange terms (lower 
right) or direct term only (upper right). 



It clearly appears that the excitation of the collective mode is due to the anti- 
symmetric T = state, while the strength arising from the symmetric T = 1 state 
lies within the particle hole continuum. 

Our correlated Hartree Fock results are in good agreement with those obtained by 
Cowell and Pandharipande using correlated states and a cluster expansion truncated 
at the two body level [31] . On the other hand, the results of Ref. [31] suggest that 
the inclusion of the exchange term of Eg. ( 15.701) is essential for the excitation of 
the coherent state, while our results show that the corresponding peak survives if 
exchange contributions are neglected. This discrepancy is likely to be ascribed to 
differences in the effective interactions. 
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Summary &; Outlook 



We have carried out calculations of the charged current weak response of nuclei and 
nuclear matter in both the low momentum transfer and the impulse approximation 
regimes. The quantitative understanding of this quantity is required in many areas of 
physics, ranging from simulations of supernovae explosions and neutron star cooling 
to the analysis of neutrino oscillation experiments. 

The calculation has been performed using a many-body approach based on a 
realistic nuclear hamiltonian, including two- and three-nucleon interactions, yielding 
a good description of the properties of both the two-nucleon systems and uniform 
nuclear matter. 

At beam energies around 1 GeV, the formalism based on the target spectral 
function has been tested comparing the calculated cross sections to electron scat- 
tering data in the kinematical region corresponding to quasielastic scattering and 
A-production, relevant to many long baseline oscillation experiments. The results 
suggest that, while the quasi elastic cross section is understood at the level of ~ 10 
%, the accurate description of the resonance region will require the extension of the 
available fits of the nucleon structure functions down to ^0.4 GeV^, as proposed 
inRef. [T03] . 

In the high energy domain the differences between our results and the prediction 
of the FG model, widely used in the analysis of experimental data, mostly arise 
from the effects of short range nucleon-nucleon correlation. On the other hand, the 
response at low |q| (of the order of tens of MeV) is known to be sizably affected by 
long range correlations, giving rise to the excitation of collective modes. 

The response associated with Fermi transitions at low momentum transfer has 
been calculated from an effective interaction, derived using the formalism of corre- 
lated basis functions and the cluster expansion technique. Our work improves upon 
existing effective interaction models in that it includes the effects of many-nucleon 
forces, which become sizable, at high density. 

The energy per nucleon of both symmetric nuclear matter and pure neutron 
matter, obtained from our effective interaction model, turns out to be in fairly 
good agreement with the results of highly refined many body calculations, based on 
similar dynamical models. A comparable agreement with the results available in the 
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literature has also been found for single particle properties, e.g. the effective mass, 
and the spin susceptibility. The emerging picture suggests that our approach can 
be regarded as an effective theory that captures the relevant physics, allowing one 
to obtain reasonable estimates of a number of different quantities using low order 
perturbation theory in the Fermi gas basis. 

The responses calculated within the correlated HF approximation show that the 
inclusion of short range correlations leads to a significant quenching of the transition 
matrix elements and shifts the strength towards larger values of the energy transfer, 
uo, for all values of |q|. At low momentum transfer, long range correlations have been 
taken into account within the TD approximation. The excitation of the coherent 
state can be clearly seen in our results at |q| = 0.3 fm~^. 

The most straightforward extension of our approach is the calculation of the 
Gamow Teller response. The inclusion of finite temperature effects, needed to extract 
the observables relevant to supernova and proto-neutron star physics, also appears 
to be doable, in the low temperature region T <^ m^. 

As a final remark, it has to be pointed out that, while our approach can and 
should be further developed, the possible improvements only pertain the structure 
of the effective interaction and the inclusion of perturbative corrections, and do not 
involve going to higher order in the cluster expansion. 

Although the contribution of clusters involving more than two nucleons is known 
to be, in general, non negligible, effective theories are in fact designed to provide 
lowest order results reasonably accounting for the available data. 

In this context, however, the rationale underlying the choice of the two-nucleon 
cluster approximation, employed in our work, deserves a comment. 

While being likely to be reasonable for use in nuclear matter limit, our procedure 
does not preserve the normalization of the many-body wave function. We plan to 
investigate quantitavely this issue using alternative truncation schemes, designed to 
guarantee the correct normalization at finite order in the cluster expansion |1081I109] . 

The most obvious improvement is the inclusion in Wgff of the non static com- 
ponents of NN potential, which are known to be needed to reproduce scattering 
data. 

On the other hand, inclusion of higher order terms in the perturbative expansion, 
which is expected to be rapidly convergent, is necessary to take into account more 
complex mechanisms, that play a role in determining several properties of many- 
body systems as, for example, the effective mass at the Fermi surface [56] . 
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Appendix A 

Properties of the operators Of- 



In this Appendix, we discuss the properties of the six operators defined in Eq. fll.l6p . 
as well as some useful properties of the Pauli matrices. 



A.l Pauli matrices 

In the standard representation, in which is chosen to be diagonal, the threee Pauli 
matrices are given by (we specialize here to the spin matrices a*: analog properties 
obviously hold for the isospin matrices r*) 



The Pauli matrices satisfy 



a'a^ = 6ij + ieijk(T^ , (A.2) 
eijkcr^a^ = 2ia^ , (A. 3) 



that can be put in the form 

[a\a^] = 2te,,ka'' , (A.4) 
{a\a^} = 26,j , (A.5) 

where e^fc is the totally antisymmetric tensor and i,j,k = 1,2,3. The first properties 
shows that the Pauli matrices are the generators of an SU (2) algebra. 
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A. 2 Projection operators 

Let now <ti and (72 be the vectors of Pauli matrices for particle 1 and 2, respectively 
(i.e. (Ti = {al,al,al}). From the properties (]A.2|) -( lA73l) . it follows that 

(O-I • CT2f = 3 - 2{(Ti ■ <T2) . (A.6) 

As (cTi ■ (T2) is a scalar quantity, we can interpret the above equation as an algebraic 
one, with solutions (cri ■ 0-2) = —3 and (cri ■ (T2) = 1. They correspond to the 
states of total spin S = (spin singlet channel) and S = 1 (spin triplet channel), 
respectively. It is thus useful introducing the operators P2S+1 (and the analog II2T+1 
for the isospin states), defined as 

P,..o,^Pi = ' - '7 • , (A.7) 
P„.,.A = l±(f^, (A.8) 

which project onto states of definte total spin or 1, respectively: 

P2s+i\S') = 5ss'\S') , (A.9) 



The projection operators satisfy to 



p2 

-^25+1 


— P2S+I 


P1 + P3 


= 1, 




= 0, 



where 1 is the two-dimensional identity matrix. 

A. 3 Spin and isospin exchange operators 

Consider the two-nucleon spin states (or the analog isospin states) 

|oo) = -L (I ID _ I , 



= lii) 
111) = ITT) 



10) = -^(|Ti) + UT) , 



(A.IO) 
(A.ll) 
(A.12) 
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A. 4 - The tensor operator S12 



where |0 0) = IS* = Ms = 0) etc., and the inverse relations 

I TT) = 111) , 

ITi) = ^(|io) + |oo)) , 

lit) = ^(|io)-|oo)) , 
I ii) = |i -1) • 

From properties ( lA.9p . and from 

(i^3-Pi)ITT) = ITT) , (A-Pi)lii) = lii), 
(P3 - Pi) I Ti) = lit) , (P3-Pi)UT) = ITi), 

it follows that P„ = — Pi is the spin-exchange operator, satisfying 

P,\S Ms) = {-f^'\S Ms) . (A.13) 

A similar exchange operator can be defined for isospin, = lis — IIi, with 

PrlTMr) = (-)^+^|TMt) . (A. 14) 

Combining the above results we find 

P,, = P,P, = 1 (^1 + (cTi . 0-2)) (1 + (ti ■ T2)) , (A.15) 

with 

P,r\SMs,TMT) = H'-^^'IS Ms,T Mt) . (A.16) 



A. 4 The tensor operator S12 

The tensor operator 5*12 is defined as 



3 

5'i2 = ^ (cTi ■ r) (cr2 ■ r) - (cti ■ erg) , (A.17) 



where r is the relative coordinate of particels 1 and 2 while r = |r|. 
Making use of Eq. flA.2|) . it can be shown that 

S'i2(cri ■ 0-2) = (cTi ■ (T2)Si2 = Si2 • (A.18) 

As we saw, (cti ■ (T2) = 1 on triplet states, while (cti ■ (T2) = —3 on singlet states. 
The above equation thus implies that the tensor operator only acts on triplet states 
and 

[5i2,P3] = 0. (A. 19) 
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Moreover, 



^1^2 = 6 - 2^12 + 2(cri ■ 0-2) 



The tensor operator is a function of r satisfying 



12 



CTi (0-2 ■ r ) + ( CTi ■ r 1 



2— (cTi ■ r) (0-2 ■ r) 



12 



60 
— ^'-'12 



(A.20) 

(A.21) 
(A.22) 



For any function u{r), Eg. flA.2ip implies 

du r 



yu) ■ (vsu) = -r-- ysi2) = 

dr r 



Moreover 



[VS 



12 



12 



36 



[^12, ySu)] = -,t{SXT 



36 



[S12, (VSu)] V = ^ (L ■ S) 



(A.23) 

(A.24) 
(A.25) 
(A.26) 



where S = (cri + (T2) /2 and L = rxp = — z(rx V) is the orbital angular momen- 
tum operator of the relative motion. 

From Equation (lA.22p . we can calculate 



and 



[S12, V^Sia] = 



ysu) [Su, V] = - ys,2y 



(A.27) 
(A.28) 



A. 5 Operator algebra 



Equations flA.6p . flA.lSP and flA.20p show that the six operators 
close an algebra, i.e. they satisfy 

k 
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(A.29) 



(A.30) 



A.6 - Matrix elements 



The coefficients K^, are easily obtained by calculating 



0^0' = O'O^ = O' =^ K^^ = i^fi = 6f 

0^0^ = 30^ - 20^ =^ K^^ = 3(5f - 2S^ , 
0^0' = 0^0^ = 0^^ = Kl^ = 6^ , 
0^0^ = 30^ - 20^ =^ K^^ = Kl^ = 5l- I5l , 

Ate = ^62 = 3^5^ - 251 , 

k 





0^0'^ = 


0^0^ 


= 0^ 




-- o'o^ = 


= 30^- 


- 20^ 




0^0' = 


30^ - 


20' - 




-- 0^0' = 


= 30^- 


-20" 




0^0^ = 


0^0^ 


= 0^ 




0^0^ = 


0^0' 


= 0' 




-60^- 


-60'- 


f40^ 




0^0'' = 


O^'O^ 


= 0' 


0^0' = 


-- o'o^ = 


= 30^- 


-20^ 


0^0^ 


= 60^ + 


20' - 


- 20^ 



Kl, = 35f - 2SI , 
. K',, = Kl, = 361 - 251 , 
' Kl, = Kl, = 51 , 

-"■36 ~ -"-63 ~ "6 ' 

> K^u = 95f - Q5l - Q5l + 4^4^ , 

> K',, = Kl, = 5l , 

> ^46 = ^64 = 355' - 256^ , 

= Q5\ + 251 - 2^5^ , 



0^0^ = 0^0^ = 60^ + 20" - 20^ =^ = K^^ = 65^ + 25l - 25^ , 



0^0^ = 180^ - 120^ + 60^ - AO" - 60^ + 40^ 



K^e = 185'l - 125^ + 6^3^ - 45^ - 65^ + 45, 



A.6 Matrix elements 

Finally, we report a number of expectation values of operators involving Pauli ma- 
trices, in two-nucleon states of definite total spin and isospion, IS* Ms, T Mt)- 

(-P25'+in2T'+l) = 5sS'^TT' , (A. 31) 

{P2S'+lIl2T'+lPar) = H'^^ 5 SS' ^TT' , (A.32) 

J2^S'l{Sl2P2S'+lT^2T'+l) = 5s'i5tT'J2^1Ms\Si2\IMs)=0 , i^-^^) 

SMs Ms 

J2Ss'l{Sl2P2S'+lli2T'+lPar) = 0. (A.34) 

SMs 
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A ~ Properties of the operators O" 



A. 7 More matrix elements 



The exphcit expressions for the matrices entering Eg. (12 .601) . defined by 



(A.35) 



where \Xfi) denotes the two-nucleon spin-isospin state, can be easily obtained from 
the above properties of the six operators O"-^. 

We find 



A' 



A' 



A' 



A' 



A' 



A' 



( 1 

1 
1 
1 

/ 1 

1 

-1 

V-1 

/ 1 

-1 
1 

/ 1 

-1 
-1 
1 

/ 1 

1 

-1 
-1 

/ 1 

-1 
-1 

V 1 

/ 1 




V 



1 
1 
1 
1 

1 
1 

-1 
-1 

-1 
1 

-1 
1 

-1 
1 
1 
-1 

1 
1 

-1 
-1 

-1 
1 
1 

-1 



1 






1 
1 
1 
1 

-1 
-1 
1 
1 

1 

-1 
1 
-1 

-1 
1 
1 

-1 

-1 
-1 
1 
1 

-1 
1 
1 

-1 




1 





1 \ 

1 
1 

1/ 

-1 \ 

-1 
1 

1 } 
-1 \ 

1 

-1 

1/ 

1 \ 

-1 

-1 

1/ 
-1 \ 

-1 
1 

1 / 
1 \ 

-1 
-1 

1/ 

\ 




1/ 



(A.36) 



(A.37) 



(A.38) 



(A.39) 



(3 cos^ B-\)=J^{Z cos^ ^ - 1) , (A.40) 



(3cos2^-l) = A'^(3cos2^-l) , (A.41) 



(A.42) 
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B' 



B' 



B' 



B' 



( 1 


2 

V 
/ 1 

2 


V 

/ 1 

2 
2 
4 

/ 1 

-1 


V 

1 

-1 

2 
-2 



v 





1 


2 

2 
1 



2 
1 
4 
2 

-1 
1 



-1 
1 

-2 
2 



2 

1 




1 
2 

2 
4 
1 
2 



1 

-1 

2 
-2 
1 

-1 



0\ 

2 


1/ 

\ 

2 

1 / 

4\ 

2 
2 

1 / 

0\ 


-1 

1/ 
-2\ 

2 
-1 

1/ 



(3 cos^ ^ - 1) 



(3cos2^ - 1) 



(A.43) 



(A.44) 



(A.45) 



(A.46) 



(A.47) 



where Q is the angle between r and the z axis. 



A. 8 Change of representation 

In this Section we discuss the different representation for the operators of the "fe" 
algebra. A generic operator x can be written as 



X = ^X^ijO^ = Xc+Xr{TrT2)+X^{cri-a2)+X„r{o-VO-2){ri-T2)+XtSi2+XtTSl2{Ti-T2) 

p=l 



in the basis of operators (]A.29p . or as 



(A.48) 



X — [xto + SsiXtrSu] -P25+in2r+i 



(A.49) 



TS 



in the "TS-representation" . 
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The transformation matrix is given by 



/ 1 -3 -3 
1 1 -3 



V 1 



-3 
1 



9 \ / Xe \ 
1 / \ Xar / 



XtT 



( a^oo \ 
V / 



xm 
xn 



or 



= + (4T - ?,)xr + {AS - 3)x^ + {AS - 3)(4T - 3)x<, 
XtT = xt + {AT - 3)xtT . 



The inverse transformation is given by 



1 

16 



3 
1 

-3 
-1 

1 3 
-1 1 



9\ 

3 
3 

1/ 



/ Xqo \ 

XlO 
Xoi 

\ ^11 / 



/ Xc \ 

X(j 



XtO 
Xtl 



Xt 

XtT 



(A.50) 



(A.51) 



(A.52) 



(A.53) 



(A.54) 
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Appendix B 

Energy at two-body cluster level 



The energy per particle at two-body cluster level can be written (see Eqs. (12.311) and 
(E33D) 

{/\E)2 = J2 I f/l2, [ ^1 + fu] 1 + /l2t;i2/l2 - Jl) , (B.l) 



with 



2 

Kj 



U = -T^V' , h + h = --V^ - -^V^ , (B.2) 
2m m 4m 



where V acts on the relative coordinate r, while V/j acts on the center of mass 
coordinate R, defined as 

r = ri-r2 , R=^(ri + r2) (B.3) 

respectively. 

Including only the static part of the interaction, i.e. neglecting the spin-orbit 
components, both the correlation function /12 and the two-nucleon potential V12 are 
written as 

6 6 
p=l p=l 

with the six operator listed in Eg. (II. 161) . whose properties are discussed in 
Appendix Rl. 

The FG two-nucleon state is given by 



) = - e'^^^-'''+^^-'''^ \S Ms,TMt) 



V 



}_ gi(k.r+K.R) \SMs,TMt) , (B.5) 
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where 

|ki|, \kj\ < pf 
k=^(k,-k,) , K = ki + kj, (B.6) 

and \S Ms,T Mt) denotes a state of total spin (isospin) S (T) and spin (isospin) 
projection Ms (Mt). 

We will discuss the potential and kinetic energy term separately. 



B.l Potential energy 

Consider the operator 

and the decomposition of /12 in the TS-representation (see Eq. flA.49p ) 



/12 = E [/■ 

ST 



ST + ^SlftTS: 



12 



P2S+1^2T+1 



(B.7) 



(B.8) 



In the above equation, -P2S+1 ^2T+i ^Pi^ isospin projection operators, 
whose properties are given in Appendix El By writing the corresponding decompo- 
sition for W12 and f 12 and calculating 



W12 



T0t^T0 + ^Sl\ 



/ri + 8/ir + 2 {fxiftT ~ frr) ^12 



TS ^ 

+ V,^ [16 {UJ,^ - f^^) + (/^, - + 12/,y ^12] } }P25+in: 



'-2T+1 ) 



we can identify 



'-'to /ro 



(/li + S/i'r) + IQv^T [friftT - ftT 

{JtiIiT " ftT^ + ^'iT {JtI ~ ^fTlftT + ^2/; 



After replacing 



Kj ij 
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(B.9) 



(B.IO) 



B.2 ~ Kinetic energy 



the potential energy contribution to {AE)2 reads 



(w) 



2V^ 



SMs TMt kikj S'T 
^S'lWtT'{'r){Si2P2S'+1^2T'+l) 



Ws'T'{r){P2S'+1^2T'+l) + 



(B.ii: 



Ws'T'{r){P2S'+1^2T'+lP<TT) + 5s'lWtT'{r){Si2P2S'+1^2T'+lPcTT) 



where P^t is the spin-isospin exchange operator defined in Appendix |A] and the 
expectation values (O) are taken over two-nucleon states of definite total spin and 
isospin I S Ms, T Mt). Using 



d^rid^r2 = I d-'r d-'R = V I d'r , 



'3„ j3 : 



the definition of the Slater function fl2.6ip . 

V 



k|<PF 



ikr 



(27r)3 



N 



\k\<pp ^ 



(B.12) 



(B.13) 



and the results of Appendix El we finally obtain 
1 1 



2V^ 1/2 



\/^(25 + l)(2T+l) f d^rwsAr) [l - (-l)^^^^^^^^^)] 

ST 



(B.14) 

where u denotes the degeneracy of the momentum eigenstates. In the case of sym- 
metric nuclear matter [u = 4) 

^(w^) = J d^r^^[woo{r) + 9wii{r)]a_{pFr) + 

+ [3woi{r) + 3wio{r)]a+{pFr)]'j , (B.15) 

where p = N/V is the density and 

a±{x) = l±f{x) . (B.16) 



B.2 Kinetic energy 



Let us now discuss the kinetic contribution to the energy, given by 



fi2, [ti +t2, /12] 



2m 



/12, [ v^ fu] 



(B.17) 



We consider spin- zero and spin-one channels separately. 
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B ~ Energy at two-body cluster level 



Spin-zero channels In these channels, the relevant part of the correlation func- 
tion is given by 

/i2 = $^/To(r)Pin2T+i . (B.18) 

T 

Making use of the results of Appendix |Al as well as of 



we find 



fro, V^/to 



fi2, [ v^ 







/to, (V/to)V 



TO) 



(B.19) 



[/to -PLn2r+i; [V^, /to] -Pin2r' 

TT' 

[/to, [v^, /to] -Pi^n27._^^n2'r,+i 

TT' 

5^ [/to, (VVto) + 2(V/to)v] AH^T+i 

T 

2$^[/to,(V/to)v] p^n 

T 

-2 5^(v/To)'Pin 



'■2T+1 



^2T+1 • 



(B.20) 



Finally, 



2m 



/12, [v^/l2]l =-5^(v/to)' An 



2T+1 • 



(B.21^ 



Spin-one channels In these channels, the correlation function is given by 



P TT 



B.22) 



fi2 = Y.[fTiir) + ftT{r)S 

T 

Relying once more on Appendix \^ we calculate 

Yl [V', (/t'I + /iT'5i2)P3n2T'+l] = E {[^'' Z^'l] + /tT'5i2]}P3n2T'+l 

= {(^'/^'l) + 2(V/iT')V + (VVtT'^12) + 2(V/tT'^12)v}P3n2T.+i 

T' 

= J2 {(^'/^'i) + 2(V/tr')V + (VVtT')'^i2 + (V'^i2)/iT' 

T' 

+ 2(V/iT')(V^l2) + 2^i2(V/iT')V + 2ftT'{VSu)V^P3ll2T'+l ■ (B.23) 
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B.2 ~ Kinetic energy 



Hence, the commutator in Eq. (]B.17p can be rewritten as 



/i2, [v^/l2] 



with 



and 



We find 
and 



[(/t1 + ftTSi2)P3 ^2T+l^ {■ ■ ■} P3 ^2T' 

TT' 

|^/t1 + ftTSu, {■ ■ •} 

T 

(-^T ^ + Ft^^P3^2T+1 > 



-P3 ng^r+i 



.(1) 



/ti,{...} 



ftTSu, {■ ■ ■} 



{(VVt'i) + 2(V/iT')V + (VVirO-^is + (V25i2)/iT' 
2{VftT'){VSu) + 25i2(V/ir')V + 2/iT'(V5i2)v} 



.(1) 



-2(VM)'-2(V/ri)(V/iT)5i2 , 



7(2) 
r 



/iT^i2,2(V/n)Vj + [/tT^i2,25i2(V/Ti)V 

+ [/iT5i2,2/Tl(V5i2)V = 
-2(V/Tl)(V/iT)5i2 - 2(V/tr)'522 + 
12 J V 



-2 (V/n)(V/tT)5i2-2(V/ir)'(8-25: 



12; 



2/. 



^ (L.S) + ^(8-5- 



12 



Collecting all pieces togheter, we find for the spin-one channels 



2m 



/12, [V^/l2]l = -VUV/Tl)' + (V/Tl)(V/iTM + 
J J Z — / I 



+ (Vftrf (8-28,2) +/, 



(B.24) 



(B.25) 



(B.26) 
(B.27) 



(B.28) 



^(L-S) + -^ (8-^12) 



-P3 ^2T+l 



TS 



'SI 



2(V/T5)(V/iT)5l2 + 
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with 



ttT 
tbT 



H^ftrfSf, + f,r^ (L ■ S) 



(8-S- 



12) 



-^25+1-^2^ 



+1 



J2 {tTs{r) + 6si [ttT{r)Su + hT{r)(L ■ S)] }p2s+i^2T+i , (B.29) 



TS 



1 ^2 



m 
1 

m 
1 



■(V/: 



TO) 



48 



{VfTlY + 8(V/tT)' + ^ //t 



6 



m 

1 36 



2(V/ti)(V/,t) - 2(V/ir)^ - - U 



m r 



2 .^tr 



B.3 Final expression 

We can rewrite 



with 



where 



iAE), = Y,{^J\Wu\^J-J^) 



i<j 



12 



m 



/i2, [ v^ 



+ /l2^12/l2 



T5 



51 



WtT{r)Si2 + W,T{r){l^-^) 



(B.30) 



-P25+in2T+l ) 



Wti 



1 



48 



+ 



-%i (/ti + S/tr) + IKt (/ri/tT - ftr) 



tT 



1 

m L 



2(V/Ti)(V/iT)-2(V/ir) 



,2 •/tr 



+ 



+ 2Vrpi(^fnpif^rp ftT^^'^tT{jTl ^friftT^^'^fi 



w, 



bT 



1 36 

,.9 



m r 



2 •'ir • 



100 



B.3 - Final expression 



Making use of the expression for the expectation values given in Appendix [XJ 
we finally obtain (compare to Eq. (lB.15|) ) 



N 



P_ 
32 

+ 



Wooir) + 9Wuir) 



a-{pFr) + 



31^01 (r) +31^10 (r) 



iB.31] 
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Appendix C 

Euler-Lagrange equations for the 
correlation functions 



C.l Spin singlet channels: uncoupled equations 

In the spin-zero channels, the energy per partcicle of SNM, evaluated at two-body 
cluster level, reads (compare to Eqs. flB.15|) and flB.211) ) 



N 



^ (2T + 1) 47r y r^dr 
/to' /to 

where aTs{x) = 1 ~ {~)'^~^^i^{x) and 



m 



aToiPpr) 
aToippr) 



const / dr F 
'o 



/to' fro 



ifrof + ^ wto/to 



with 



(C.l) 



(C.2) 



(C.3) 



The corresponding Euler-Lagrange (EL) equations for the unknown functions /to 
are given by 

(C.4) 



From 



dF 



TO 



2 m vto fro <Pto 



TO 
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dF 

Wo 
d dF 



2 /to 'P^o ! 



2/: 



TO 'PtO + 4 /to '/'to 'i'^ro 5 



we obtain 



Introducing 
we can put Eq.( ]C.6|) in the form 

fi'ro " 



dr (9/^0 

/to 4>to + 2 /to 0TO - ^ ^To fro 0to = 

fl'TO = /to 0rO ; 



^TO 

6^0 



+ m t^To fl'TO = 



(C.5) 

(C.6) 
(C.7) 

(C.8) 



Now we introduce a Lagrange multiplier, in order to fulfill the requirement (see 

Eqs.(123S])-(I23H])) 

gTo\r=d = (i>'To\r=d ■ (C-9) 



The resulting equation is Eq.(4) of Ref. |107j 



9to 



V^TO 

0ro 



+ m {vto + A) 5-^ = 



to be integrated with the boundary conditions 



5^701^=0 
9T0\r=d 



(C.IO) 



(C.ll) 
(C.12) 



C.2 Spin triplet channels: coupled equations 



In the spin-one channels, the contribution to the energy is given by (see Eqs. llB.lSP 
and dEnD) 



(A^)^ 

N 



yfrif + 8(V/,^)^ + 5 f^^ 



+ 



const / dr F 

'0 



/ti' At! /ti5 ftT 
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(C.13) 



C.2 - Spin triplet channels: coupled equations 



where 



/ti5 ftr'^ /ti5 ftT 



(/ti) V^i + 8 (/^) V^i + ^ A't0ti + 



+ m 



VTl{f^l + 8f^T)+16VtT{fTJtT-ftT) ■ (C.14) 



In this case we have two coupled EL equations 







d dF 


dF 


dr a/^^ 




d dF 


dF 


dr a/j'j. 


df„ 



. 

Carrying out the derivativees as in the spin-zero channels and defining 

9X1 = /ti0ti , QtT = VSftT(pTl , 

we find 

9ti - (j^ + ^ ^Ti) gri - mVSvtTgtT = 

^ + m (vti - 2vtT) + ^ gtT- mVSvtTgTi = . 
Finally, inclusion of the Lagrange multipliers needed to guarantee 

gTl\r=di ^ 0Tllr=di ' 
gtT\r=d2 ~ 't'Tl\r=d2 ' 

with, in general, di ^ d2, leads to (compare to Eq.(5) of Ref. |107j ) 

g'Tl- + (^Tl + Ai) gTl-m[\/%VtT + >^2)gtT = ^ 



g'lT 



(C.15) 



(C.16) 



(C.17) 



(C.18) 
(C.19) 



g'lr 



1^ + m {vti - 2vtT + Ai) + ^ gtT - m {VSvtT + A2) gri = 



with the boundary conditions 



gTi\r=o 

9Tl\r=di 



and 



gtT \r=0 
gtT\r=d2 



= 0, 
= . 



(C.20) 

(C.21) 
(C.22) 



(C.23) 
(C.24) 
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